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Abstract 

Over a long period of exploration, the successful observation of quantized version of anomalous 
Hall effect (AHE) in thin film of magnetically-doped topological insulator completed a quantum 
Hall trio—quantum Hall effect (QHE), quantum spin Hall effect (QSHE), and quantum anomalous 
Hall effect (QAHE). On the theoretical front, it was understood that intrinsic AHE is related to 
Berry curvature and U(l) gauge field in momentum space. This understanding established con¬ 
nection between the QAHE and the topological properties of electronic structures characterized by 
the Chern number. With the time reversal symmetry broken by magnetization, a QAHE system 
carries dissipationless charge current at edges, similar to the QHE where an external magnetic field 
is necessary. The QAHE and corresponding Chern insulators are also closely related to other topo¬ 
logical electronic states, such as topological insulators and topological semimetals, which have been 
extensively studied recently and have been known to exist in various compounds. First-principles 
electronic structure calculations play important roles not only for the understanding of fundamen¬ 
tal physics in this field, but also towards the prediction and realization of realistic compounds. In 
this article, a theoretical review on the Berry phase mechanism and related topological electronic 
states in terms of various topological invariants will be given with focus on the QAHE and Chern 
insulators. We will introduce the Wilson loop method and the band inversion mechanism for the 
selection and design of topological materials, and discuss the predictive power of first-principles cal¬ 
culations. Finally, remaining issues, challenges and possible applications for future investigations 
in the field will be addressed. 
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I. INTRODUCTION 


I.A. Hall Effect and Anomalous Hall Effect (AHE) 

Edwin H. Hall discovered, in 1879, that when a conductor carrying longitudinal current 
was placed in a vertical magnetic held, the carriers would be pressed towards the transverse 
side of the conductor, which led to observed transverse voltage. This is called Hall effect 
(HE)^, and it was a remarkable discovery, although it was difficult to understand at that 
time since the electron was not to be discovered until 18 years later. We now know that 
the HE is due to the Lorentz force experienced by the moving electrons in the magnetic 
held, which is balanced by a transverse voltage for a steady current in the longitudinal 
direction. The Hall resistivity p^y under perpendicular external magnetic held H (along the 
z direction) can be written as p^y = RqH, where Rq is the Hall coefficient, which can be 
related to the carrier density n as Rq = —^ (in the free electron gas approximation). The 
HE is a fundamental phenomenon in condensed matter physics, and it has been widely used 
as as experimental tool to identify the type of carrier and to measure the carrier density or 
the strength of magnetic helds. 

In 1880, Edwin H. Hall further found that this “pressing electricity effect” in ferromag¬ 
netic (EM) conductors was larger than in non-magnetic (NM) conductors. This enhanced 
Hall effect was then called as the anomalous Hall effect (AHE)^, in order to distinguish it 
from the ordinary HE. Later experiments on Ee, Co and Ni^ suggested that the AHE was 
related to the sample magnetization M (along z), and an empirical relation for the total 
Hall effect in EM conductors was established as^’^ 

Pxy — RqH + RsM, (1) 

where the second term is the anomalous Hall resistivity, and its coefficient Rg is material- 
dependent (in contrast to Rq, which depends only on carrier density n). 

As schematically shown in Eig. 1, if we plot the Hall resistivity p^y versus external 
magnetic field H, we will generally expect a straight line (which crosses the origin of the 
coordinates) for NM conductors (see Eig. 1(a)); however, non-linear behavior will be ex¬ 
pected for EM conductors, wherein p^y increases sharply at low field and crosses over into 
a linear region under high field (see Eig. 1(b)). The initial sharp increase of p^y is due to 
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FIG. 1: (Color online) Typical behaviors of the Hall and anomalous Hall effects. The Hall 
resistivity p^y is plotted versus external magnetic field H. (a) The Hall effect; (b) The 
anomalous Hall effect; and (c) The hysteresis loop measured from the anomalous Hall 

effect. 

the saturation of magnetization of the sample under external field. After the saturation of 
magnetization, p^y will depend on H linearly (for high field region), which is dominated by 
the ordinary Hall contribution. Therefore, the slope of the linear part under high field gives 
us Rq. If we extrapolate this linear part to the zero field limit (H=0), it will not go through 
the origin of coordinates, and its intercept on y-axis gives us RgM as can be learned from 
Eq. (1). This is a surprising result, which immediately implies two important facts: (1) 
A kind of Hall effect can be observed even in the absence of external magnetic field (i.e., 
no Lorentz force); (2) The anomalous Hall resistivity is sensitive to magnetic moment M, 
and it has been suggested that this property would be useful for the detection of magnetiza¬ 
tion of conducting carriers, particularly for the cases of weak itinerant magnetism, such as 
surface/interface magnetization®, dilute magnetic semiconductors^’®, etc. Experimentally, 
measurements are usually done in magnetization loops by scanning the magnetic field from 
positive to negative, and a hysteresis loop in the p^y vs. H relationship is to be expected. 
This is very similar to the familiar hysteresis loop observed in the M vs. H curve, as usually 
found in FM conductors. The empirical relation Eq. (1) is very simple and widely used. 
Unfortunately, its correctness is not well justified. As will be pointed out in a latter part 
of this article, the anomalous Hall resistivity may have a very complicated form, which is 
usually non-linear in M. 

Although the HE and AHE look quite similar phenomenologically, their underlying 
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physics are completely different. The HE is due to the Lorentz force’s influence on the 
moving electrons under magnetic field; however the AHE exists in the absence of external 
magnetic field, where there is no orbital effect of moving electrons. In other words, why do 
the electrons move towards the transverse direction in the absence of Lorentz force? The 
mechanism of AHE has been an enigmatic problem since its discovery, a problem that lasted 
almost a century. This problem involves concepts deeply related to topology and geometry 
that have been formulated only in recent years®“^^ after the Berry phase was recognized in 
1984 ^®. 

Karplus and Luttinger provided a crucial step in unraveling this problem as early as 
1954^^. They showed that moving electrons under external electric field can acquire an 
“anomalous velocity”, which is perpendicular to the electric field and contributes to the 
transverse motion of electrons and therefore to the AHE. This “anomalous velocity” comes 
from the occupied electronic states in EM conductors with spin-orbit coupling (SOC). They 
suggested that the mechanism leads to an anomalous Hall resistivity proportional to the 
square of the longitudinal resistivity: p^y oc p^.^. Since this contribution depends only on 
the electronic band structures of perfect periodic crystal and is completely independent 
of scattering from impurities or defects, it is called intrinsic AHE, and it was not widely 
accepted until the concept of Berry phase was well established. 

Eor a long time, two other extrinsic contributions had been considered to be the dominant 
mechanisms that give rise to the AHE. Smit^®’^® suggested that there always exist defects 
or impurities in real materials, which will scatter the moving electrons. In the presence 
of SOC and ferromagnetism, this scattering is asymmetric and should lead to unbalanced 
transverse motion of electrons. This is called skew scattering, and Smit argued that it is 
the main source of the AHE. This mechanism predicted p^y ^ Pxx, in contrast with the 
intrinsic contribution. Berger^^, on the other hand, argued that electrons should experience 
a difference in electric field when approaching and leaving an impurity, and this leads to 
another asymmetric scattering process called side jump, which also contributes to the AHE. 
This mechanism marvelously predicts Pxy ~ p^.^, the same as the intrinsic contribution. 

The debates about the origin of AHE lasted for a long time, and no conclusion could 
be drawn unambiguously. From the experimental point of view, defects or impurities in 
samples are unavoidable and are usually complicated with rich varieties. The contributions 
from both intrinsic and extrinsic mechanisms should generally coexist^^. Experimentalists 
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have tried hard to distinguish them^^ The controversy also arises because of a lack of 
quantitative calculations that could be compared with experiments. 

Still, the discovery of the quantum Hall effect in 1980s^® and the later studies on the 
geometric phase^^ and the topological properties of quantum Hall states^*“^°, promoted the 
fruitful study of AHE significantly. Around the early years of this century, the intrinsic 
AHE mechanism proposed by Karplus and Luttinger was completely reformulated in the 
language of Berry phase and topology in electronic structures®^^^. It was recognized for 
the first time that the so-called “anomalous velocity” originates from the Berry curvature 
of occupied eigen wave functions, which can be understood as effective magnetic field in 
momentum space. This effective magnetic field modifies the equation of motion of electrons 
and leads to the intrinsic AHE®^^^. Erom this understanding, we now know the intrinsic 
AHE can be evaluated quantitatively from the band structure calculations, thanks to the 
rapid development in the field of first-principles calculations. Detailed calculations for var¬ 
ious compounds portray the contributions from intrinsic AHE in a way that convincingly 
confirms the existing experimental results^^’^^’^^’^^’^^’^^, establishing the dominant role of the 
intrinsic contribution to the AHE. In the meantime, the extrinsic AHE was also analyzed 
carefullyand it was understood that its contribution dominates for the clean limit case 
(rather than the dirty limit), where the deviation of distribution of electronic states from the 
equilibrium distribution is significant. In this article, we will not discuss the extrinsic AHE. 
Eor readers who want to learn more details, please refer to the review article by Nagaosa 
et al.^^ Introduction of the Berry phase mechanism for the understanding of intrinsic AHE 
was a big step forward in the field, and it is also the fundamental base for understanding 
the quantum anomalous Hall effect (QAHE), which is the main subject of this article. 


I.B. Berry Phase Mechanism for Intrinsic AHE 


In quantum mechanics, the Berry phase is the quantal phase acquired by the adiabatic 
evolution of wave function associated with the adiabatic change of the Hamiltonian in a 

parameter space 7 ?., with 'R. = (7^i,7^2,.being the set of parameters (a vector). 

Let \n(7t)) be the n-th eigenstate of the Hamiltonian The overlap of two wave 
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functions infinitesimally separated by AT?- in the 7?.-space can be evaluated as 

(n(7?.)|n(7?. + AT?.)) = 1 + AT?.(n(T?.)|V 7 ?,|n(T?.)) = exp[—iAT?. • An(T?.)], (2) 

where A„(T?-) = i(n(T?-)|V 7 ;.|n(T?.)) is called the Berry connection. Here A^(T?-) is an 
important quantity, because it can be viewed as a vector potential and its curl 0^(T?.) = 
Vt?, X A„ (T?.), called the Berry curvature, gives an effective magnetic field in the parameter 
space T?. (as will be addressed below). The Berry phase 7 „ can then be defined as the 
integral of the Berry connection along a closed loop C in the parameter space, or according 
to Stoke’s theorem, equivalently as the integral of Berry curvature on the surface S enclosed 
by the adiabatic loop C = dS (i.e., the surface with the loop C as boundary): 

7„ = / AJT?.) ■dlZ= [ ■ dS, (3) 

Jc Js 

where the second equality suggests that the Berry phase can also be regarded as the effective 
magnetic flux passing through a surface S. 

Although the concept of the Berry phase has broad applications in physics, its relevance 
to the band structure in solids has been recognized only in limited situations, such as the 
quantum Hall effect under a strong magnetic field^® and the calculation of electronic polariza¬ 
tion in ferroelectrics^^’^^. Here we will show that the Berry phase concept is also important 
for understanding intrinsic AHE. In this case, we treat the crystal momentum (k) space as 
the parameter space. We consider a crystalline solid with discrete translational symmetry. 
The eigen equation of the system is given as H{r)'i/;n-k{r) = £nkV'nk(i’); where H{r) is the 
Hamiltonian, and V'rik(r) are the eigen energy and the eigen wave function, and n is the 
band index. Because of the translational symmetry, the eigen wave function ipnki'^) of fho 
system is k-dependent, where k is the momentum defined in the first Brillouin zone (BZ) of 
momentum space. According to Bloch’s theorem, the eigen wave function should take the 
form of Bloch state V’nk(r) = e**^“’u„k(r), with u„k(r) being the periodic part of the wave 
function. Then the eigen equation of the system can be recast as i7k(r)unk(i’) = £nk'ifnk(r), 
where iTk = is the k-dependent Hamiltonian. Following the above discussion, 

we can now define the Berry connection and Berry curvature in the parameter (momentum 
k) space as 

■^n(f^) ^(^nk|^k|^nk) 

^n(f^) ^k X A^(k) ^(VkUjjkl X IVk'fin.k)' 
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( 4 ) 

( 5 ) 


These two quantities are crucial in understanding the Berry phase mechanism for intrinsic 
AHE. Before we can go further, we have to clarify several important properties of the Berry 
connection and Berry curvature. 

• Berry connection is gauge dependent: As we have learned from the textbooks of solid 
state physics, there exists an arbitrary phase factor for the eigen wave function |urik) 
that is not uniquely determined by the eigen equation of the system. Under a U(l) 
gauge transformation, the eigen wave function [unk) is transformed into |u^k) 

l«nk> = (6) 

where 0n(k) is a real and smooth scalar function of k. It is easy to see that |u^k) i^ 
the eigen wave function of the system for the same eigen state, i.e., i^kl^^nk) = ^nkl^^nk) 
is satisfied. However, the corresponding Berry connection will be changed by such a 
gauge transformation, becoming 

A;(k) = i{<k|VU<k> = A„(k) - Vk.;.„(k). (7) 

If we regard the scalar function 0„(k) as a kind of scalar potential in the k-space, 
this form of transformation is the same as that of a vector potential of magnetic 
field in real space, and therefore the Berry connection A„(k) can be viewed as an 
effective vector potential in momentum space. The gauge dependence of the Berry 
connection suggests that it is not physically observable. However, it becomes physical 
after integrating around a closed path (i.e., the Berry phase 7 „ defined in Eq. (3)). 
This is because the integration of the second term of Eq. (7) around a closed path will 
only contribute an integer multiple of 2tt, and the Berry phase is therefore invariant 
modulo 27r. 

• Berry curvature is gauge invariant: This conclusion can be drawn directly from the 
factor that Vk x Vk0n(k) = 0; therefore r2„(k) = Vk x A(j(k) = Vk x A„(k) is 
unchanged under the U(l) gauge transformation of Eq. (6). This vector form of Berry 
curvature suggests that it can be viewed as an effective magnetic field in momentum 
space. It is a gauge invariant local manifestation of the geometric properties of the wave 
function in the parameter (k) space, and has been proven to be an important physical 
ingredient for the understanding of a variety of electronic properties^^’^^’^^’^^’^^. 


• Symmetry consideration: Here we consider two important symmetries in solid state 
physics, namely time reversal symmetry (TRS) and inversion symmetry (IS). Following 
the above discussion, it is easy to learn that the Berry curvature has the following 
symmetry properties: 


= 12n(-k) for inversion symmetry, 

ll„(k) = -0„(-k) for time reversal symmetry, 


( 8 ) 

( 9 ) 


which suggests that 0„(k) = 0 for a system with both TRS and IS. They are strong 
symmetry constraints, which means that in order to study the possible physical effects 
related to the Berry curvature, a system with either broken TRS or broken IS is 
generally required. For example, the intrinsic AHE is observed in a ferromagnetic 
(FM) system with SOC, where the TRS is broken. For a system with TRS only, in 
general r2„(k) ^ 0, however its integration over the whole Brillouin zone should be 
zero as implied by Eq. (9). In such a case, the Berry curvature may take effect only 
through k-selective or band n sensitive probes. 

• The ehoiee of periodie gauge: For the convenience of studies on solid states with 
discrete translational symmetry, we usually use the periodic gauge in practice. That 
is, we require the Bloch wave function to be periodic in momentum space, and use 
the condition t/^k(r) = ^k+G(r) (where G is a translational vector for the reciprocal 
lattice) to partly fix the phases of wave functions. This periodic gauge is used in most 
of the existing first-principles electronic structure calculation codes. However, we have 
to note that it is a very weak gauge-fixing condition, and abundant gauge degrees of 
freedom still remain. 

• Gauge field: Spreading out each component of the vector form discussed above (Eq. 
(5)), the Berry curvature can be written explicitly as an anti-symmetric second-rank 
tensor 






( 10 ) 


9 









where is the Levi-Civita antisymmetric tensor, and /r, v, ^ can be simply treated 
as the direction indices (i.e., x,y,z) for our purposes. For example, for many of our 
following discussions on the transport properties in a two dimensional (2D) system, 
we will use the equality. 


( 11 ) 



Writing it in such a way also helps us understand that Berry curvature is simply the 
U(l) gauge field in the language of gauge theory. In the presence of gauge 

freedom, in order to preserve the gauge invariance of Lagrangian under the gauge 
transformation, we need to use the gauge covariant derivative instead of the original 
derivative. For example, for the position operator = idk^, its gauge covariant form 
is written as x^ = idk^ — A^{k), where A^^k) is the vector potential (i.e., the Berry 
connection). In such a case, the gauge covariant position operators x^ and Xy are 
no longer commutable, and their commutation relation is given as [x^^x^] = 
which gives the gauge field = dk^Ay - dk^A^. 

Keeping in mind all these properties discussed above for the Berry connection and Berry 
curvature, we are now ready to understand the intrinsic AHE from the viewpoint of gauge 
field and geometric Berry phase. For simplicity, let us first consider a 2D system with 
Hamiltonian Ho{x, y) defined in the xy-plane. Under an external electric field Ey, the whole 
Hamiltonian of the system becomes H{x,y) = HQ{x,y) + eEyy. By evaluating the velocity 
operator as x = — |[x, iif] = — |[x, i^o] ~ ^Eyy], we find that the presence of electric field 
in the system leads to an additional velocity term, ^Eyy] — —^ExyEy = —^ilzEy, 

which is along the x direction (being transverse to the external electric field Ey). This result 
comes from the non-vanishing commutation relation between the gauge covariant position 
operators. This additional velocity term is nothing but the anomalous velocity term proposed 
by Karplus and Luttinger^^, and is the origin of the intrinsic anomalous Hall effect. 

Now let us generalize the above discussion to three dimensional (3D) space and consid¬ 
ering both the electric and magnetic external fields, we find that the equation of motion of 
electrons, in the presence of gauge field, should be recast as^® 



( 12 ) 

( 13 ) 
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where E and B are the external electric and magnetic fields, respectively. It is apparent 
that the group velocity Eq. (12) is different from the textbook one by the inclusion of an 
additional contribution (i.e., the second term on the right-hand side). This term, k x On, 
is now known as the anomalous velocity term, and it is proportional to the k-space Berry 
curvature r2n(k). This equation of motion can also be obtained from the semiclassical theory 
of Bloch electron dynamics. To do this, we have to consider a narrow wave-packet made 
out of the superposition of the Bloch state of a band. Here we will not discuss the detailed 
derivation of this approach; readers please refer to the review article by Xiao et al. for more 
details^^. 

For a ferromagnetic system in the presence of SOC, the time reversal symmetry is broken; 
therefore the non-vanishing Berry curvature contributes to the anomalous velocity term. 
Consider a system under external electric field Ey ^ 0 and without magnetic field (i.e., 
B = 0). The current density can be then obtained from the equation of motion shown 
above as 


= -ev^ = -ey^ 

n 

-EL 


d^k 


BZ ( 271 ") 

d^k 


h 


E 


BZ (27r)' 

d^k 


BZ 


(2vr) 


3/„(k)(-kxl2Jk)), 
/„(k)(ExO„(k)), 
fn{k)EyQ.n,z{^), 


(14) 


where fn{k) is the Fermi distribution function. The dc Hall conductivity is then given as 

— r 




^xy — 

H/qi 


e 

h 




BZ 

27ih 


-fn{k)nn,z{k) 


E r /' (/n(k)f!„,.(k)) 

^ J —TV j —TV 


dkxdky 


dkz 


(15) 


(16) 


which is simply the Brillouin zone (BZ) integral of the Berry curvature weighted by the 
occupation factor /„(k) of each state. For the convenience of our later discussions, in the 
last part of the equation, we have explicitly written the 3D integral as a 2D integral in the 
^x^ 2 /-plane, supplemented with a line integral along kz- 

On the other hand, we know from textbooks that the DC Hall conductivity can be also 
derived from the Kubo formula as 




m^n 


'BZ 


d^k 

(^ 


(/n(k) - um X Im 


I I ^^mk) {'^mk \ \ "^nk) 

(^nk ^mk)^ 


(17) 
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which is basically equivalent to Eq. (15), by expanding each term of the vector form of the 
Berry curvature explicitly with some algebra. However, we have to emphasize that the form 
of Eq. (15) for Hall conductivity is not only compact, but also fundamentally important 
for the underlying physics. It demonstrates the Berry phase mechanism of the intrinsic 
AHE, and it is the key equation for our following discussions on the geometric meaning and 
topological nature of QAHE. Eor example, for a two dimensional system, the BZ integral 
of Berry curvature for a fully occupied band must give rise to an integer multiple of 27r. 
Therefore, the Hall conductivity a^y for a 2D insulator can be finite and be an integer 
multiple of jh. This is the underlying physics for the quantization of the AHE. 


In addition to the understanding of the Berry phase mechanism for the intrinsic AHE 
discussed above, another important step forward in this field is the quantitative and accu¬ 
rate evaluation of the anomalous Hall conductivity (AHC) for realistic materials, yielding 
results that can be compared with experimental observations. Due to the rapid increase 
of computational power and the development of first-principles calculation methods, such 
quantitative comparisons (though still difficult) become possible and contribute greatly to 
the development of the field. In a series of papers reporting quantitative first-principles 
calculations for SrRuOs^^, Fe^^ and CuCr 2 Se 4 _a;Bra;^^, it was demonstrated that the calcu¬ 
lated AHC can be reasonably compared with experiments, suggesting the importance of the 
intrinsic AHE. In more recent years, accurate calculations of AHC have been achieved by 
using the Wannier function interpolation scheme^^’^®. The presence of SOC and the breaking 
of TRS are crucially important for the intrinsic AHE in those systems. 


In summary of this part, we have discussed the Berry phase mechanism for intrinsic AHE. 
Conceptually, this is an important increment in our understanding because the intrinsic 
AHE is now directly linked to the geometric and topological properties of the Bloch states 
in momentum space. This understanding has also led to a great number of research projects 
on the exotic properties of electronic systems with SOC, such as QAHE and the topological 
insulators. 
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II. QUANTUM HALL FAMILY AND RELATED TOPOLOGICAL ELECTRONIC 
STATES 


One of the main subjects of condensed matter physics is the study and classification of dif¬ 
ferent phases and various phase transitions of materials that possess rich varieties of physical 
properties. Landau developed a general symmetry-breaking theory to understand the phases 
and phase transitions of materials. He pointed out that different phases really correspond 
to different symmetries of compounds. An ordered phase, in general, can be described by a 
local order parameter, and the symmetry of the system changes as a material changes from 
one phase to another. Landau’s symmetry-breaking theory is very successful in explaining a 
great many kinds of phases (or states) in materials, but it does not work for the topological 
quantum states. A topological quantum state cannot be described by a local order parame¬ 
ter, and it can change from one state to another without any symmetry-breaking. Topology, 
a word mostly used in mathematics, is now used to describe and classify the electronic struc¬ 
tures of materials. “Topological electronic state” means an electronic state which carries 
certain topological properties (in momentum space usually), such as the states in topological 
insulators^^^^^, topological semimetals®^^^^, topological superconductors^®’®°’^^^™, etc. One 
of the most important characteristics of topology is the robustness against local deforma¬ 
tions, or in the language of physics, the insensitivity to environmental perturbations, which 
makes topological electronic states promising for future applications. To characterize the 
order of topological states, new parameters, namely the topological invariants, are needed. 
In this section, we will review some of the topological electronic states and their quan¬ 
tum physics, including the integer quantum Hall (IQH) state, the quantum anomalous Hall 
(QAH) state, the quantum spin Hall (QSH) state, and the topological semimetal (TSM) 
state. We will use the concept of Berry phase in momentum space to discuss the topological 
nature of topological states and the related topological invariants. 

II.A. Chern Number, Chern Insulator, and Quantum AHE 

In 1980, Klaus von Klitzing discovered^® that the Hall conductance, viewed as a function 
of strength of the magnetic field applied normal to the two-dimensional electron gas plane, 
at very low temperature, was quantized and exhibited a staircase sequence of wide plateaus. 
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The values of Hall conductance were integer multiples of a fundamental constant of nature: 
e^/h = 1/(25812.807572 fl), with totally unanticipated precision, and independent of the 
geometry and microscopic details of the experiment. This is an important and fundamental 
effect in condensed matter physics, and it is called the integer quantum Hall effect (IQHE) 
— the quantum version of the Hall effect. 

The IQHE is now well understood in terms of single particle orbitals of an electron in 
a magnetic filed (i.e., the Landau levels), and the phenomenon of “exact quantization” has 
been shown to be a manifestation of gauge invariance. The robustness and the remarkable 
precision of Hall quantization can also be understood from the topological nature of the 
electronic states of 2D electron gas under magnetic field. The integer number, originally 
known as the TKNN number^® in the Hall conductance derived from the Kubo formula, 
is now characterized as a topological invariant called “Chern number”. This topological 
understanding of the IQHE is a remarkable leap of progress, opening up the field of topo¬ 
logical electronic states in condensed matter physics. IQHE is therefore regarded as the first 
example of topologically non-trivial electronic states to be identified and understood. 

This conceptual breakthrough with regard to the topological nature of IQHE, though 
important, was not widely generalized, and the relationship of IQHE to the rich variety of 
condensed materials was not revealed for a long time. This is because the IQHE is observed 
only in a very particular system, i.e. 2D electron gas under strong external magnetic field, 
and the formation of Landau levels (usually at very low temperature) is required. Under 
such extreme conditions, the material’s details and its electronic band structure become 
irrelevant to the physics. In this sense, the lattice model proposed by Haldane in 1988®^ is 
very stimulating for the study of topological electronic states. He proposed that a spinless 
fermion model on a periodic 2D honeycomb lattice without net magnetic flux can in principle 
support a similar IQHE. Although Haldane’s model is very abstract and unrealistic (at least 
at that time), his result suggested that a kind of quantized Hall effect can exist even in 
the absence of magnetic field and the corresponding Landau levels. In other words, certain 
materials, other than the 2D electron gas under magnetic field, can have topologically non¬ 
trivial electronic band structures of their own, which can be characterized by a non-zero 
Chern number. Such materials were called Chern insulators later. In this way, the concept 
of topological electronic state was generalized and was connected to the electronic band 
structure of materials. 


14 


The progress in the study of IQHE was very inspiring in the 1980s; on the other hand, 
not much of it was related to the field of AHE. Around that time, the development of the 
AHE field was rather independent, with almost no intersection with the study of IQHE. But 
the understanding of Berry phase mechanism for intrinsic AHE®^^^, reformulated around 
the early years of the 21st century, changed the situation significantly. This understanding 
established the connection between the AHE and the topological electronic states. By this 
connection, it is now understood^°’^^’®° that intrinsic AHE can have a quantum version — 
the quantum anomalous Hall effect (QAHE), which is in principle similar to IQHE but 
without external magnetic field and the corresponding Landau levels. It turns out that this 
is exactly the same effect discussed by Haldane for a Chern insulator with a non-zero Chern 
number. Currently, the fields of IQHE and AHE are closely combined, and the realization of 
a Chern insulator becomes possible by finding the QAHE in suitable magnetic compounds 
with strong SOC. 

The IQHE and the QAHE have their own characteristics; however, their underlying 
physics, in terms of the topological properties of their electronic structures, are basically 
the same. They are all related to the Berry connection. Berry curvature and Berry phase in 
momentum space, as discussed in the previous section. In this part, we will not discuss the 
particular details associated with each effect, but rather we will concentrate on their common 
features, namely the topological property and topological invariant (Chern number). To 
begin with, we will discuss several important concepts: 


• Hall Conductivity as Integral of Berry Curvature: Considering a 2D insulating system 
with broken TRS, the Hall conductivity of the system at low enough temperature can 
be written as 
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where the second equality follows because of the existence of an energy gap and /n(k) = 
1 (=0) for the occupied (unoccupied) state (see Eqs. (14)-(16) for a 3D version). This 
form of Hall conductivity can be derived either from the Berry phase formalism or 
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FIG. 2: Under the periodic boundary condition for the Bloch states, the 2D Brillouin zone 

is equivalent to a torus. 

from the Kubo formula, as discussed in the previous section. It is also important to 
note that this formula can be equally applied to the anomalous Hall conductivity of a 
2D FM insulator without external field and to the Hall conductivity of a 2D electron 
gas under magnetic field. For the former it is straightforward; for the latter case, 
however, special considerations have to be taken. First, we assume that the external 
magnetic field is strong enough for the system to form the Landau levels. Second, 
the system must be insulating with the chemical potential located within the gap 
(i.e., between two neighboring Landau levels). Third, since the presence of external 
magnetic field will break the translational symmetry, we have to use the concept of 
magnetic translational symmetry to recover the periodicity of the system. As a result, 
magnetic Bloch states must be used for the evaluation of Berry curvature, and the 
magnetic Brillouin zone is used for the integration, correspondingly. 

• 2D Brillouin Zone as a Torus: By adopting the periodic gauge, the wave function 
'i/ink(r) is periodic in momentum space. Under such a periodic boundary condition, 
the 2D BZ becomes a torus as schematically shown in Fig. 2. Therefore, the integral 
of the Berry curvature over the 2D BZ can be recast as a vector integral over the 
surface of the torus as 


Here we have used an effective Berry curvature 0(k) defined over the torus surface, 
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which satisfies the condition 0(k) • n = r 22 (k), where n is the unit vector along the 
normal of the surface S. 

• Quantization and Chern Number: The torus is a closed manifold without boundary, 
so according to Chern’s theorem®^ in differential geometry, the vector integral of Berry 
curvature over a torus surface must be an integer multiple of 27r, 

J 0(k) • dS = 27rZ, (20) 

where Z is an integer number called Chern number. Having this property, it is now 
easy to see that the Hall conductivity for a 2D insulator must be quantized as 

=-r 27rZ = ——Z. (21) 

^ 27r/j /j ^ ^ 

• Absence of Smooth Gauge: Based on Stokes’ theorem, the surface integral of Berry 

curvature (Eq. (20)) can also be evaluated as an loop integral of the Berry connection 
A(k) along the boundary of the BZ, i.e., J^Sl(k) • dS — ■ dk. Since the BZ 

is a torus, which has no boundary, the integral must be vanishing if A(k) is smoothly 
defined in the whole BZ. Therefore, a non-zero Chern number indicates that we cannot 
choose a smooth gauge transformation such that A(k) is continuous and single valued 
over the whole BZ. Thus, non-zero Chern number can be viewed as obstruction to 
continously define the phase of the occupied wave functions on a 2D BZ, which is a 
torus.®^“*^. 

For a system with multiple bands, the Berry curvature should be understood as the 
summation of contributions coming from all occupied bands. Having the properties discussed 
above, we can now define a Chern insulator as a 2D insulator whose electronic structure gives 
a non-zero Chern number. The above discussions about the Chern number and the topology 
of electronic states in momentum space may be still too abstract, in the following we will 
give a more explicit explanation, in which the none zero Chern number manifests itself as the 
winding number of the Wannier center evolution for the effective ID systems with constant 
ky. Here we consider the simplest case: a 2D lattice with only one single occupied band 
(the band index n can therefore be neglected). Although it is impossible to choose a single 
smooth gauge over the whole BZ for a Chern insulator according to our above discussions. 
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it is possible to choose a special gauge that is smooth and periodic along one direction (say 
kx), but not necessarily along the other (say . Thus, we can do the integration in Eq. 

(19) explicitly for the k^ direction with fixed ky and obtain (we set the lattice parameter 
a=l) 
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FIG. 3: The evolution of 0{ky) (or ID Wannier center) on cylinder surface for system with 

(a) Z=0, (b) Z=1 and (c) Z=2. 

Here, 0{ky) = dk^A^ is an angle (i.e.. Berry phase) calculated from the ID integration 
of Ax{kx, ky) along the k^ axis (a closed loop) for each fixed ky. We can then plot 0{ky) as 
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a function of ky. As shown in Fig. 3, over a cylinder surface (in the cylinder coordinates), 
ky is along the longitudinal direction, and the azimuth is the angle 0{ky) for each fixed ky. 
Moving from ky=—'K to ky='K, we can see the difference between the trivial insulator {Z=0) 
and the Chern insulator (Z 7 ^ 0). The winding number of 0{ky) over the cylinder surface is 
zero for the former (Fig. 3(a)), and non-zero for the latter (Fig. 3(b), (c)). In this way, we 
relate the Chern number to some kind of winding number defined by the eigen functions, 
and see that the Chern insulator has a “twisted” energy band. 

It is also interesting to note that the fc^-dependent Berry phase calculated as 

PniK) = ^ h) ( 23 ) 

can be related to the Wannier center of a ID system (where we have recovered the band 
index n dependence). To see this, we consider a generic 2D system and denote the Wannier 
function in cell R associated with band n, given in terms of the Bloch states, as 

"" bAh / , (24) 

(27r)^ Jbz 

with 

|^nk) = 5]e'‘^-^|nR), (25) 

R 

where k = {kx,ky) is considered, and the periodic part of the Bloch function is defined as 

«nk(r) = e“*‘"''’?/^„k(r). ( 26 ) 


Then the matrix elements of the position operator r between Wannier functions take the 
form 

{nR\r\mO) = Vkl^mk), (27) 

(27r)^ Jbz 

from which we get the center of the Wannier function in the 0 -th cell as 


r = (n 0 |f|n 0 ) = ^ 77 ^ / d^k(M„k|Vkl^nk)- 
(27r)^ Jbz 


(28) 


Now let us consider the k^ and ky directions separately, and define the ID Wannier 


function for each fixed-as 
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Then the ID Wannier center (along x direction) as function of ky can be defined as the 
average value of position operator x as 

Mnk) = Pn{ky), (30) 

from which we see that Pn{ky) is nothing but the ID Wannier center of band n. The 
Chern number can therefore be understood as the winding number of ID Wannier center 
when it evolves as a function of ky (see Fig. 3). We can also see that when the Chern 
number is not zero, we cannot choose a smooth gauge transformation such that A(k) is 
continuous and single valued over the whole BZ. Although the above discussion is only for 
the simplest systems with single occupied band, generalization to the multi-band situation is 
quite straightforward. For a general band insulator with multiple occupied bands, the Berry 
connection A„m(k) = f (unkl Vklumk) contains band index n, m and becomes nonabelian 
Then the Berry curvature that defines the Chern number is obtained from the trace of the 
Berry connection A(k). 

The Hall conductivity of a Chern insulator with a non-zero Chern number must be 

2 

quantized as an integer multiple of y. Different Chern numbers give different states whose 
Hall conductivities are different, but their local symmetry can be the same. Therefore, to 
distinguish the states, we cannot use local order parameters (in the language of Landau’s 
symmetry breaking theory); instead we need a topological invariant, the Chern number, as a 
global order parameter of the system. By rewriting Hall conductivity in terms of the Berry 
curvature and Berry phase, we can now unify QHE and QAHE. Readers will recall that the 
QAHE is nothing but the quantum version of the AHE realized in a Chern insulator without 
the presence of external magnetic field. 

Similar to QHE, where the TKNN number is related to the number of edge states in a real 
2D sample with boundary^®, the Chern number can also be physically related to the number 
of edge states for a 2D Chern insulator^^. The existence of edge states is a direct result 
of the topological property of the bulk electronic structure, and is due to the phenomenon 
discussed in the literature as the bulk-boundary-correspondence®®’®^. In this case, due to the 
broken TRS, the edge state must be chiral (i.e, the electrons of the edge state can move only 
in one direction surrounding the sample boundary, either left- or right-handed, as shown in 
Fig. 4(b)). As discussed with regard to IQHE, the charge transport of the edge state is in 
principle dissipationless, and back-scattering is absent due to the lack of an edge state with 


x{ky) = {n, Rx = 0, ky\x\n, Rx = 0, ky) = 


(27r) j_, 


dkx {Un\i IV kx 
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FIG. 4: Schematic illustrations of the edge states in (a) QHE; (b) QAHE, and (c) QSHE. 


opposite velocity®®’®^. 

Although the topological properties of Chern insulators and the related QAHE are fun¬ 
damentally the same as that discussed with regard to IQHE, they are conceptually broadly 
generalized to a wide field and to a rich variety of materials. This generalization is an impor¬ 
tant step forward, providing the building blocks for subsequent discussions of many possible 
topological electronic states. Erom the application point of view, the Chern insulator (or the 
related QAHE) is also important because the quantized Hall conductivity can be realized in 
the absence of magnetic field, greatly simplifying measurement conditions. 

II.B. Z2 Invariant, Topological Insulator, and Quantum Spin Hall Effect 
(QSHE) 

Although the Chern insulator (or the related QAHE) is the simplest topological electronic 
state, its realization occurred rather later and was much stimulated by the rapid develop¬ 
ment in the field of topological insulators—another interesting topological electronic state 
protected by TRS. Considering symmetry, it is easy to prove that the Chern number in Eq. 
(20) must be vanishing for an insulator with TRS. However, this does not mean that the 
electronic state carries no topological property in this case. Kane and Mele^® introduced a 
new topological invariant, the Z 2 number, to classify an insulating system with TRS. They 
proposed that a time reversal invariant insulator can be further classified as a trivial in¬ 
sulator with Z 2 = 0 or a non-trivial topological insulator with Z 2 = 1, which is a typical 
example of a symmetry protected topological state®°. A Z 2 topological insulator in 2D is 
also called a quantum spin Hall insulator (QSHI) because it can support the quantum spin 
Hall effect (QSHE)^®“^^, which shares certain features with the IQHE and QAHE and can 
be understood from the viewpoint of “band-twisting” or winding number of Berry phase 
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in momentum space. A QSHI is different from a trivial insulator in the sense that it has 
gapped insulating states in the bulk but gapless states on the edge due to its topological 
property. It is also distinguished from the Chern insulator in the sense that QSHI has 
even number of edge states, composed of pairs of counter-propagating chiral edge states, 


as shown in Fig. 4(c)). The QSHE has been explicitly discussed for graphene lattice^® and 
HgTe quantum well structures^°“^^. In this part, however, we will take a simple task and 


discuss the invariant and the topological state with TRS from the viewpoint used above 
for the Chern insulator. 

For a 2D insulator with TRS, the total number of occupied electronic states must be 
even. Suppose we can divide them into two subspaces, I and II , which are related by TRS 
and smoothly defined on the whole BZ. Evaluating the Chern number for each subspace 
independently, if one has Chern number Z, the other one must have Chern number -Z (due 
to the TRS) and the total Chern number of the whole system is always zero. It seems 
that we can use Z as the topological index to classify the band insulators with TRS. While 
unfortunately, the smooth partition of the occupied states into two subspaces with one being 
the time reversal of another is only possible with extra good quantum numbers, such as 
where the Chern number obtained within the spin up subspace can be used to describe the 
topology of the system and is called “spin Chern number”. But generically such a smooth 
partition of the occupied states can only be made for half of the BZ (not the whole BZ) 
and we need at least two patches (A and B) to fully cover the whole BZ. The winding 
number of the U(2N) gauge transformation matrix tAsik) at half of the boundary between 
two patches can be used to define a new classification of the band insulators with TRS. As 
further proved by Eu and Kane, only the even and odd feature of the above winding number 
is unchanged under the U(2N) gauge transformation and all the band insulators with TRS 
can be classified into topological trivial and non-trivial classes, which is called Z 2 invariance 
accordingly. Fu and Kane further derived the expression for Z 2 invariance in terms of both 
Berry connection and Berry curvature as®^ 



(31) 


where integral of Berry curvature is performed over the half BZ (i.e., HZ/2), and the integral 
of Berry connection is performed along the boundary of the half BZ (i.e., d{BZ/2)). To 
evaluate Z 2 invariance using the above equation, one needs to find a smooth gauge for the 
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wave functions on half of the BZ, which is very difficult for the band structure calculations for 
realistic materials, and as a result, this formula is rarely used in practice. On the other hand, 
the Z 2 invariance is a gauge invariant quantity, one should be able to compute it without 
any gauge fixing condition. For that purpose, some of the authors of the present paper 
developed an alternative expression for the invariance called Wilson loop method, which 
will be introduced in detail in section 3.2. Here we just sketch its main idea briefly. Similar 
with the Wilson loop method for Chern number calculation, the purpose of the Wilson loop 
method is to calculate the “Wannier center” of each band for the effective ID insulators with 
fixed ky and determine the Z 2 invariance by looking at their evolution with ky. With the 
presence of TRS, a generic system contains 2N occupied bands. As mentioned previously, 
now the Berry connection becomes 2N*2N matrix. The loop integral of such U(2N) Berry 
connection gives a 2N*2N unitary matrix D{ky). The U(l) part of this matrix contributes 
to the Chern number as discussed before and the Z 2 invariance can be obtained from the 
remaining SU(2) part by taking the phases of its eigenvalues 9n{ky), where n denotes the 
band index. As proved in Ref. 92, the TRS only guarantees the double degeneracy of 6n{ky) 
at two time reversal invariant loops ky = Q and ky = tt. The Z 2 invariance of the system is 
determined by whether the On{ky)s switch parters or not as they evolves from ky = 0 to tt. 

To be explicit, let us again consider the simplest example of 2D insulator, with only 
two occupied bands, V'/.k and which are related by the TRS (say, for example, 

where T is the time reversal operator). The two states form a Kramers 
pair and must be degenerate at the time-reversal-invariant momentum (TRIM) of the BZ 
(i.e, with kx or ky = 0,7r). We can evaluate 0(ky) for all occupied bands (as we have done 
above for the Chern insulator). As mentioned previously, the two Wannier centers must be 
degenerate ai ky = 0 and tt, which leads to three situations in general (as shown in Fig. 
5(a)-(c)). First, if the evolution path of the 0{ky) doesn’t enclose the half BZ, this is the 
trivial situation with Z 2 invariance 1^=0 (Fig. 5(a)). Second, if the evolution path of the 
0{ky) enclose the half BZ once, this is the non-trivial topological insulator with Z 2 index 
u=l, where the crossing of two 0{ky) curves at ky = 0 and tt is protected by the TRS (Fig. 
5(b)). Third, if the evolution path of the 0{ky) enclose the half BZ twice, the two 0{ky) 
curves must cross at some ky other than the TRIM. Such crossings are not protected by the 
TRS and are removed by small perturbations, which drive the crossings into anti-crossings. 
As the result, the system becomes equivalent to the trivial case with Z 2 index u =0 (Fig. 
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FIG. 5: Evolution of Wannier centers (Pn{ky)) for systems with (a) ^2=0, (b) ^2=1, (c) 
\Z\=2 but ^2=0. The Wannier centers wind the cylinder surface even (a, c) and odd (b) 

times respectively. Figure from Ref. 92 


5(c)). 

The edge states of 2D Z 2 topological insulators must appear in pairs due to the existence of 
TRS, and each pair of edge states is composed of two counter-propagating chiral edge states, 
which are related to each other by the TRS. In this case, if one state at the edge has velocity 
V, we can always find another state with opposite velocity —v. However, back-scattering 
between the two states is again forbidden because the two states must have opposite spins. 
Existence of such pair of edge states leads to the observable where the spin 

Hall conductivity is quantized in units of ^ (if Sz is conserved). According to our above 
discussions of the Z 2 topology, if there exist multiple pairs of edge states, two pairs can in 
principle couple together and open up a gap (because scattering between the pairs is not 
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forbidden). As a result, the Z 2 number is physically related to the number of pairs of edge 
states mod 2 , in contrast to the case of Chern insulators. 

Using the 2D topological invariants (such as Z or Z 2 ) as building blocks, a 3D crystal 
can in principle be characterized by a triplet of 2D topological numbers for the three crystal 
orientations, respectively. In such a way, we can extend the topological classes from 2D to 
3D. There are in general two situations by doing this: (1) one is the trivial extension leading 
the “weak” 3D topological class, which can be viewed as a simple stacking of 2D topological 
insulating layers along a certain crystal orientation; ( 2 ) another is a non-trivial extension 
leading to a “strong” topological class, which is a new and stable topological state in 3D. 
For example, by a simple stacking of 2 D Chern insulating layers along the direction, if the 
band dispersion along is weak enough, each 2D layer with fixed k^ will keep its original 
topological property (characterized by a non-zero Chern number). As a result, a 3 D “weak” 
Chern insulator can be obtained. In a similar way, a 3 D “weak” topological insulator can 
be obtained too. However, it is important to note that, if we concentrate on insulators only, 
there will be no “strong” Chern insulator in 3D; on the other hand, a new 3D state with 
TRS, called “strong” topological insulator, can be obtained®^“®^ by extending the Z 2 number 
from 2 D to 3 D. (As will be discussed in the next section, the missing “strong” Chern class 
in 3D has to be a metal—the topological semimetal). 

For a 3 D insulator with TRS, there will be 8 TRIM points in the BZ (different from 2 D BZ 
which has only 4 TRIM points). For any 4 TRIM points that lie in the same plane, we can 
define a 2D Z^ number as discussed above. In total, there will be 16 invariant configurations 
and can be distinguished by four independent Z^ indices, UQ^UxUyPz), as discussed in Ref. 93. 
Here i^o is the total Z 2 number, and the i'x,y,z are the 2 D Z 2 numbers along the x,y,z 
directions, respectively. The three weak Z 2 indices i'x,y,z can be evaluated for the kx — tt, 
ky = TT, and kz = tt planes, as a convention, respectively. The total Z 2 index i/q, however, is 
a strong index, and it has to be evaluated globally by considering the 3D structure. That is, 
we have to consider the change of 2D Z 2 number between two parallel planes (perpendicular 
to certain reciprocal lattice vector). For example, we can choose the kx=0 and tt planes, or 
ky=0 and tt planes, or kz=0 and tt planes. If the Z 2 numbers of two planes are different, 
we have 1^0 = !, otherwise, if they are the same, we have 1/0 = 0 . A realistic example will 
be given in Section 3. For a 3D strong topological insulator with strong topological index 
i/q 7 ^ 0, on its boundary (i.e, surface), we will expect the Dirac cone type surface states, 
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which are protected by the TRS^®’™. 


II.C. Magnetic Monopole, Weyl node and Topological Semimetal 

As mentioned above, if we attempt to consider the possible “strong” topological state 
characterized by Chern numbers in 3D, we will hnd that it has to be a metal, and this leads 
to a very interesting new topological state of quantum matters — topological metals or 
semimetals. We will see that this is a much more general state, and the 2D Chern insulating 
state can be regarded as a special cut of the 3D topological semimetal state along a certain 
plane in momentum space. From another point of view, we can also raise a question like 
the following. Since we know from the above discussions that insulators can be further 
classified as topologically trivial and non-trivial insulators, can we do the same classification 
for metals? If possible, what will be the topological invariant for the proper description of 
the topological metallic state? This is a very interesting question, and we will show in this 
part that it is related to the magnetic monopoles^^’^® and Weyl nodes in momentum space. 

Berry phase is called quantum geometric phase since it has a very intuitive geometric 
picture. It is proportional to the solid angle subtended at the magnetic monopole by the 
adiabatic loop of Hamiltonian in momentum space. The magnetic monopole is the source 
or drain of the gauge field 0,i(k). Since 

n^k) - /m ^ _ , (32) 

it diverges at the point k = ko where Em0^o)=En0^o) (see Eqs. (16) and (17)). This means 
that the magnetic monopole is formed by energy level crossing, and a two-energy-level system 
is the simplest case. Consider the generic form of the Hamiltonian for any two-level system, 

i7 = do(k)J 2 x 2 + cl(k)-(T, (33) 

where 12x2 is the identity matrix, d is a 3D vector depending on momentum k, and a are 
Pauli matrices. There are two eigen states xp± with eigen energies E± = do(k) i \/d • d. 
The term do(k) is just a shift of zero energy level, and can be neglected. At the energy 
degeneracy or level crossing point (k = ko), E^=E_ (i.e, d(ko)=0) is required. Obviously, 
ko is not necessarily on the path of a adiabatic loop change of the Hamiltonian, which means 
that the band gap between E^ and E_ can be well preserved if the Fermi level is away from 
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the band-crossing point. Around the neighborhood of the ko point, d(k) can be expanded 
as d(k) d(ko) -h (k — ko) • Vd(k). Taking the zero point of the parameter space as ko = 0, 
we have d(k)=(k • V)d(k). Now for the simplest case, where d(k) = ±k, we can have 
i7(k) = ±k • a, and the corresponding Berry curvature can be obtained as 


0±(k) 



(34) 


Obviously, such a magnetic field distribution in momentum space is similar to the electric 
field distribution of a point charge in real space, and can be understood as magnetic field 
around a “magnetic charge”-magnetic monopole^h In other words, the divergence of mag¬ 
netic field r 2 (k) is no longer zero (i.e., Vk ■ fi(k) ^ 0) but is related to the magnetic charge 
±1 at the source or drain. The degeneracy or level crossing point kg is the place where the 
magnetic monopole with strength | is located. 

Joshua Zak®^ pointed out that such Berry phase and Berry curvature can also exist in 
periodic systems where the eigenstates are Bloch wave functions and the parameter space 
is the crystal momenta k which can vary in closed loops, such as Brillouin zone or Fermi 
surface. Then, the adiabatic evolution loop forms a compact manifold that has no boundary. 
For example, the Brillouin zone of a 2D system is a torus, and the Fermi surface of a 3D 
metal (in the simplest case) is an enclosed sphere. Obviously, Gauss’s law ensures that the 
total flux penetrating the closed surface must be quantized and is equal to the magnetic 
charge of monopoles enclosed by the surface (either torus or Fermi sphere). 

Now we can see that the magnetic monopole in momentum space (or parameter space 
in general) plays a crucial role in determining the topology of electronic band structures. 
Besides the insulating case, in the following we concentrate on a 3D metal with a well-defined 
Fermi surface in momentum space. Taking the above two-level system as an example and 
assuming that di(k) = k^, d2(k) = ky and d3(k) = k^, the magnetic monopole is located at 
the origin k^ = ky = kz = 0 and the low-energy Hamiltonian can be written as 


iJ(k) = ±k ■ a 



kx 'iky kz 


(36) 


This Hamiltonian was first proposed by Weyl®®, who found that, for massless fermions, 
the 4x4 Dirac representation is reducible, and is composed of two 2x2 (irreducible) Weyl 
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fermions with positive (+) and negative (-) chirality (and opposite magnetic charge). Similar 
to the definition of Chern number for a torus (in a 2D insulator), the total flux of the gauge 
field passing through a Fermi surface of 3D metal must be quantized as a multiple of 27r. 
We can then define the Fermi surface Chern number Cps as 


Cfs = ^ [ 0(k) • dS, (36) 

27r Jfs 

and use it as a new topological invariant to describe topological metallic states. The Fermi 
surface Chern number Cps is non-zero if a Weyl node (or a magnetic monopole) is enclosed by 
a Fermi surface, and this leads to a non-trivial topological metallic state (called Weyl metal). 
If the Fermi level happens to be exactly at the Weyl node, we will get a topological semimetal 
state (i.e, Weyl semimetal). Unfortunately, according to the “no-go theorem”for any 
lattice model, the Weyl nodes with opposite chirality have to appear in pairs (although they 
might be separated in momentum space), and the summation of Cps for aU pieces of Fermi 
surfaces must be vanishing^^^. This makes it difficult to give a proper definition for Weyl 
metal. Nevertheless, we should note that Weyl nodes (and magnetic monopoles) are stable 
topological objects, which can be well defined in 3D momentum space^°^’^°^. Two Weyl 
nodes with opposite signs may be separated in momentum space and lead to two pieces of 
Fermi surfaces, each of which has a non-zero Fermi surface Chern number Cps ^ 0. The 
two Weyl nodes may annihilate each other if and only if they overlap in k-space^^’^^“®®’^°^. 

Special attentions must be paid to the systems with both TRS and IS, where Kramers 
degeneracy exists for every momentum k. In such cases, each pair of Weyl nodes, if any, 
must overlap exactly in the k-space. In other words, the minimum effective Hamiltonian 
to describe such system must be at least 4x4 and contain two Weyl nodes (with opposite 
signs) simultaneously, as shown below. 


/f(k) 


k ■ cr 0 
0 —k • cr 


(37) 


This is called a 3D Dirac node, and is a straightforward extension of the 2D graphene 
to 3D space^°^. As we have discussed above, in such a case, a perturbative mass term can 
be introduced in principle, opening up a gap and leading to an insulating state. However, 
if we consider additional symmetries in the system, such as the crystalline symmetry, the 
mass term may again be forbidden, and this would stabilize the 3D Dirac cone and lead to 
a 3D Dirac metal or semimetal state^®’®°. Therefore, the 3D Dirac node is not as stable as 
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the Weyl node, however, it is a good starting point for us to reach the true Weyl semimetal 
state by breaking either TRS or IS in the 3D Dirac semimetals^®. 

In general, the energy level crossing may happen in many materials and at any energy 
and momentum; however, it is quite rare to have band-crossing exactly at Fermi level, par¬ 
ticularly for the cases with broken TRS or IS. Up to now, only a few materials have been 
theoretically proposed to host such Weyl semimetaU®^®^’^°® or Dirac semimetal states®®’®®’®^. 
Among them, only NasBi®®’®^ and Cd 3 As 2 ®^’®^ have been confirmed experimentally to be 
Dirac semimetals. The proposed Weyl semimetals breaking TRS, such as pyrochlore Iri- 
date^^ and HgCr 2 Se 4 ^®, have not been confirmed yet due to experimental difficulties. The 
proposals for Weyl semimetals keeping TRS but breaking IS are thought to be a way to 
overcome those difficulties. Presently the following are typical proposals: One is a super¬ 
lattice system formed by alternately stacking normal and topological insulators^®®. The 
second involves Tellurium, Selenium crystals or BiTel under pressure^®®’^®^. The third one 
is the solid solutions LaBii_ 2 :Sba,Te 3 and LuBii_a,Sba,Te 3 ^®® tuned around the topological 
transition points^®^. The fourth is TaAs-family compounds, including TaAs, TaP, NbAs and 
NbP, which are natural Weyl semimetals and each of them possesses a total of 12 pairs of 
Weyl points^®®’^®®. 

The Weyl semimetal is a new state of quantum matters. It is of particular interest here 
that the Weyl semimetal with broken TRS is closely related to the Chern insulators and 
provides a unique way to realize QAHE. Let us first consider a Weyl semimetal with a single 
Weyl node in a continuous model defined as H{k) = VF\^ cr , which as mentioned previously 
generates a monopole in Berry curvature right at the origin. The Gauss’s theorem then 
requires that the total flux flowing through out any closed surface enclosing the Weyl point 
must be ±1, i.e, equal to the chirality of the Weyl point (Eq. (36)). For a lattice system, 
there is an important “no-go theorem” indicating that Weyl points must appear in pairs 
with opposite chirality. Originally this theorem has been proved by field theory.®®’^®® Here 
we provide a much intuitive way to understand it in terms of the concept of Chern number. 
Given a Weyl semimetal in 3D lattice (assuming all trivial states are far away from the Fermi 
level), first let’s consider a 2D plane with fixed A:^, where band structures within the plane 
should be fully gapped unless the plane cut through a Weyl point exactly. The integral 
of Berry curvature over such a 2D insulating plane must be quantized and gives rise to a 
well-defined integer Chern number. Moving the plane (i.e, kz) from -vr to vr, we will then get 
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the Chern number as a function of kz, as shown in Fig. 6. Now, it is important to note that 
the Chern number as a function of kz should jump by +1 or -1 (depending on the chirality 
of the Weyl point), whenever the moving-plane goes across a Weyl point. This is because 
there exists a topological phase transition at the Weyl point, and the band gap of 2D plane 
is closed and re-opened when the plane goes across the Weyl point. This jump can also be 
understood by the following consideration. Selecting two planes at the opposite sides of a 
Weyl point, we can construct a closed manifold surrounding the Weyl point, i.e, the cube 
formed by the two parallel planes and the four side-surfaces, shown in Fig. 6 as the shaded 
area. The flux flowing through the four side-surfaces should exactly cancel each other, being 
net zero. Then, the total flux flowing through the closed manifold (i.e, the cube), which is 
now equivalent to the difference of Chern numbers for the two parallel planes, must be equal 
to the chirality of the Weyl point enclosed by the manifold (namely ±1) as discussed above. 
Having understood the jump, let’s look at the Chern number evolution as the function of kz 
from -TT to TT. In order to satisfy the periodic boundary condition in the BZ of lattice, the 
positive jump (-1-1) and the negative jump (-1) must appear in pairs, if there is any. This 
leads to the important conclusion that Weyl points with opposite chirality must appear in 
pairs. 

Finally, the total Hall conductivity of the system is given by the integral of kz as 

= (38) 

In such a way, we see the relationship between 2D Chern insulators and 3D Weyl semimetals. 
Given a 3D Weyl semimetal as defined above, we can form a 2D thin film (or quantum well 
structure) such that kz is quantized by the sample thickness, and only particular values of 
kz are allowed. If it happens that, for some particular thicknesses, the 2D Chern number 
are non-zero, we will expect quantized total Hall conductivity, and this leads to a 2D Chern 
insulator and the QAHE. 

To explicitly show the Hall conductivity of a FM Weyl semimetal, we introduce a simple 
Hamiltonian based on a cubic lattice BHZ modeT^°, 

= sm{kx)(Jx + sm{ky)ay -|- [2 -|- — cos{kx) — cos{ky) — cos{kz)](Jz- (39) 

It has two Weyl nodes at (0,0, ±fc^) with cos(k'^)=es when 0 < e* < 1. The energy dispersion 
along kz axis is shown in Fig. 7 for es=f. The Hall conductivity <Txy{kz) is calculated by 
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FIG. 6: Two Weyl nodes with opposite chirality (±1) distribut along in 3D momentum 
space. The Chern number Z{kz), evaluated for each 2D plane with fixed kz, is shown at the 
bottom as a function of kz- Jumps of Chern number are seen when the plane moves across 
a Weyl point. A cube containing one Weyl node is shadowed. 


using Eqs. (18) and (32) for each plane with fixed kz- When chemical potential fip is 0 and 
passes exactly through the Weyl nodes, (7xy{kz) is quantized to be 1 ior \kz\ < k^ and 0 for 
\kz\ > kz-, consistent with the above discussions based on the continuum model. Integrating 
<^xy{kz) over kz, the total Hall conductivity a^y is given as which depends only on 

the separation of Weyl nodes in momentum space (i.e, fc^). Shifting the chemical potential 
jjp away from 0 leads to partially occupied bands, and the corresponding Fermi surfaces 
change from two isolated spheres {0 < fip < 0.5) to two touched ones {fip = 0.5) and finally 
they merge into a peanut shape (0.5 < iip < 1.5), as shown in Fig. 7. The total Hall 
conductivity decreases with rising chemical potential (ip, because the Berry curvature from 
the upper band tends to cancel out the contribution from the lower band. Finally, when 
both bands are fully occupied the sum of Berry curvature over all bands must vanish^^^. 


31 

























Chemical potential jlp 


FIG. 7: (a) The kz band dispersion of FM Weyl semimetal, (b) The AHC dependence on 
the position of chemical potential jip- Inset: Fermi surface for different hf- 


It is important to note that a fully occupied band (without a corresponding Fermi surface) 
may still contribute to the Hall conductivity, and this should be carefully considered in 
practical calculations. In such a case, the Hall conductivity is usually evaluated by the 
volume integrals of the Berry curvature of occupied bands over the whole BZ, as is done for 
the calculation of intrinsic this is also discussed by Chen et al.^^^ and Vanderbilt 

et al.36,ii5-ii7_ 
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III. FIRST-PRINCIPLES CALCULATIONS FOR TOPOLOGICAL ELECTRONIC 


STATES 

Theoretical predictions, particularly first-principles electronic structure calculations 
based on density functional theory^^®, have played important roles in the exploration of 
topological electronic states and materials. This is not an accidental success, but rather 
due to several deep reasons: (1) To describe complicated electronic band structures of real 
materials, particularly with the spin-orbit coupling (SOC) included, first-principles calcula¬ 
tions are necessary; (2) First-principles calculations nowadays can reach the accuracy even 
up to 90% for most of the physical properties of “simple” materials (i.e., weakly correlated 
electronic materials), which makes prediction possible; (3) The topological electronic prop¬ 
erties of materials are robust, non-perturbative, and not sensitive to small error bars, if any. 
Given those advantages of first-principles calculations and its great success in this field, 
however, we have to be aware that the numerical determination of topological invariants 
(such as integer Z or Z 2 numbers) is still technically demanding, because those numbers are 
related to the phases of eigen wave functions, which are gauge dependent and randomized 
in most of the calculations. In addition to that, we also have to note that: (1) We still 
have a band-gap problem in either the local density approximation (LDA) or the general¬ 
ized gradient approximation (GGA) for the exchange-correlation potential; (2) The present 
first-principles calculations based on LDA or GGA can not treat strongly correlated systems 
properly; (3) The evaluation of Berry phase^^ and topological numbers may require some 
additional complicities, such as fine k-points meshes, gauge-fixing condition, etc. In this 
section, we will discuss some important issues and techniques related to the first principles 
studies of topological electronic states. 

III. A. Calculations of Berry Connection and Berry Curvature 

From the computational point of view, a self-consistent first-principles electronic structure 
calculation for a real material is typically performed based on the LDA or GGA for the 
exchange-correlation potential. Such a calculation will generate a set of single particle eigen 
states with wave functions Unk(i’), eigen values £ri(k), and occupation /„(k) for the ground 
state of the system. Having obtained those quantities, the remaining task is to evaluate the 
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Berry connection A„(k) = ^(MnklVkl^nk) and the Berry curvature 0„(k) = Vk x A„(k). 
The Hall conductivity can be obtained either from the summation of Berry curvature over 
the BZ (Eq. (16)) or from the Kubo formula directly (Eq. (17)). In the latter case, the 
matrix element of velocity operator v is given as 

1 d 1 d 

(t^mkl^lt^nk) ^ \ [ j 77k] l^nk) (k) ) {^mk. \ l^nk) • (^0) 

The strategy is simple, but the computational task is hard, because although the Berry 
curvature does not depend on the gauge choice, the Berry connection does. It is also true, ac¬ 
cording to our above discussions, that choosing a smooth gauge for a topologically non-trivial 
state is generally difficult or even impossible. In most calculations, the Berry connection is 
a rapidly varying function of momentum k, which makes the convergence of the calculation 
difficult. To overcome this problem, a straightforward way is to increase the number of 
k-points in momentum space. Another possible way is to use the Wannier representation, 
such as the maximally localized Wannier functions^^®. 

The first-principles calculations of Berry curvature and its integration over the BZ to 
get the intrinsic anomalous Hall conductivity (AHC) have been performed for SrRuOs and 
body-centered cubic Ee by Fang et al.^^ and Yao et al.^^, respectively. These were the 
first quantitative demonstrations of strong and rapid variation of Berry curvature over the 
BZ. The sharp peaks and valleys make the calculation of AHC challenging. Millions of k 
sampling points are required to ensure convergence. Wang et al.^^ proposed an efficient 
approach employing the Wannier interpolation technique. They first perform the normal 
ab initio self-consistent calculation by using relatively coarse fc-grid — good enough to ensure 
the total energy convergence. Then, Maximally Localized Wannier Functions (MLWF) are 
constructed for the states around the Fermi level. The number of Wannier functions should 
be carefully chosen so that the isolated group of bands around the Fermi level can be well 
reproduced, especially those bands passing through the Fermi level (because the small energy 
splitting due to spin-orbit coupling or band anti-crossing can produce sharp peaks or valleys 
of Berry curvature). Once MLWFs are obtained, the interpolation technique is applied to 
obtain the necessary quantities, such as eigen energies, eigen states, and velocity operator, on 
a much finer A:-grid. Since the number of MLWFs is typically much smaller than the number 
of basis functions used in ab initio calculations, this approach saves a lot of computation 
cost, both in time and in storage. 
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In additional to the direct methods discussed above, some advanced techniques with well- 
controlled convergency with respect to the random phases can be developed. This is because 
Berry connection is not the physical quantity that we can measure, and our ultimate goal is 
to evaluate the Berry curvature which is gauge-invariant. For example, it was proposed that 
for calculation of electrical polarization within the modern theory of polarization^^’^^ the 
change of polarization is directly proportional to the sum of Berry phase terms of occupied 
states. Since Berry phase is the integral of the Berry connection along a closed loop in the 
parameter space (see Eq. (3)), it is well defined and gauge-independent. To calculate the 
change of polarization along the direction of . 2 , a series of loops (called Wilson loops^^°’^^^) 
in momentum space are constructed. Each loop has fixed \^i_={kx, ky), but is periodic along 
kz- We can then discretize kz as ki=jGz (j=0, 1, ..., J-1), where Gz is the reciprocal lattice 
vector along kz- Taking the periodic gauge (k^,fcO), the Berry phase 

along each loop can be calculated by 


7(kj_) = Im 


j-i 


In det (k_L ) I (k_L ,fci+i)) 


j=Q 


(41) 


with band indices n and m for occupied states. By doing this, the random phase problem 
associated with eigen wave functions can be largely eliminated, and the main computational 
task is to evaluate the inner product of two neighboring wave functions along the loop. 

Eor a metallic system, the intrinsic AHC can also be calculated through this method, 
as has been implemented by Wang et al.^®. It was justified by Haldane that the non- 
quantized part of intrinsic AHC in a EM metal can be recast as a Eermi surface integration of 
the Berry curvature, so we can construct the Wilson loops over the Fermi surface rather than 
the Brillouin zone. Considering a FM metal with magnetization M along the direction, a 
2D plane (in momentum space) perpendicular to kz will cut through the Fermi surface. The 
intersection between the plane and the Fermi surface defines loops, which are the Wilson 
loops used for calculating the AHC. Following our above discussions, we can discretized 
each loop and evaluate the Berry phase along the loop efficiently by using Eq. (41). To do 
this, several things should be borne in mind: (1) The periodic boundary condition has to be 
adapted as Un,(kx,A:/) = if fhs loop does not cross the Brillouin zone boundary; (2) 

The direction of Wilson loops should be defined consistently; (3) Eor the case with multiple 
branches, a continuity condition must be used to choose each branch carefully. Finally, 
in additional to the non-quantized part of AHC, the quantized contribution to the AHC 
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(coming from the deep occupied states) can be determined by the Fermi-sea integration of 
the Berry curvature, which is gauge invariant and has no ambiguity 

This Wilson loop method is a convenient tool, and it can also be applied to investigate 
the topological invariant for systems with or without TRS as will be discussed below. 

III.B. Wilson Loop Method for Evaluation of Topological Invariants 

In real materials, the band structures are usually very complicated, and the determina¬ 
tion of topological invariants becomes essential and computationally demanding. The band 
degeneracy, either accidentally or due to symmetry, makes the numerical determination of 
phases of wave functions a tough task. Here we will present the Wilson loop method for 
determining topological indices efficiently 85,92,122,123 method is computationally easy 

and is equally applicable for Chern insulators, Z 2 topological insulators (2D and 3D), and 
the topological crystalline insulators^^^’^^^, which have attracted lots of interest recently. 

For band insulators with inversion symmetry, the indices can easily be computed as the 
product of parity eigen values for half of the occupied states (Kramers pairs have identical 
parities) at the time-reversal-invariant-momentum (TRIM) points^^®. In such a case, parity 
is a good quantum number for the TRIM points and can be computed from the eigen 
wave functions obtained from the first-principles calculations. This method has been used 
frequently and is very efficient. Unfortunately, for general cases where inversion symmetry 
is absent, the evaluation of the Z^ number becomes difficult. Three different methods have 
been considered for such a case: (i) Directly compute the Z 2 numbers from the integration 
of the Berry curvature Sl(k) over half of the Brillouin zone^^^. In order to do so, one has 
to set up a fine mesh in the k-space and calculate the corresponding quantity for each k 
point. Since the calculation involves the Berry connection A(k), one has to numerically fix 
the gauge on the half BZ, which is not easy for the realistic wave functions obtained by 
first-principles calculations, (ii) Start from an artificial system with inversion symmetry, 
and then adiabatically deform the Hamiltonian towards the realistic one. If the energy gap 
never closes at any point in the BZ during the deformation process, the realistic system 
must share the same topological nature with the initial reference system, whose Z^ number 
can easily be counted by the parity eigenvalue formula^^®. Unfortunately, making sure that 
the energy gap remains open on the whole BZ is very difficult numerically, especially in 3D. 
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(iii) Calculate the boundary (edge or surface) states. Due to the open boundary condition, 
the first-principles calculation for the boundary states is numerically heavy. The Wilson 
loop method, which we will present here, has the following advantages: firstly, it uses only 
the periodic bulk system; second, it does not require a gauge-fixing condition - thereby 
greatly simplifying the calculation; third, it can easily be applied to a general system with 
or without inversion symmetry. 

As discussed in the previous section for 2D insulators, to determine the topological num¬ 
ber (either Z or Z 2 ), the key is the angle and its ID evolution. In the single band 

model, the 0{ky) angle is given simply as the ID integration of ky) along the k^ axis. 

However, for real compounds with multiple bands, its computation requires the Wilson loop 
method®^. Once the 9 angle is computed, all the topological numbers can be determined 
easily, following our previous discussions. For the mathematic details of this method, please 
refer to our paper®^. Here we will focus on the practical steps for real calculations. 

Considering an insulator with N occupied states, after the self-consistent electronic struc¬ 
ture calculations, we get the converged charge density, the total energy, the band struc¬ 
ture, etc. Then, without loss of generality, we consider a periodic plane of {k^, ky) with 
—TT < kx < TT and —tt < ky < tt, and consider the k^ axis as a closed loop. For each fixed 
ky, we discretize the k^ line into an W mesh with mesh-index i ranging from 0 to W “ 1- 
Due to the periodic condition, the i=Nk point is equivalent to the i=0 point. The Wilson 
loop method consists of the following steps: 

(1) We need to calculate the eigen wave functions of occupied states for all mesh points. 

This can be done by continuing from the converged charge density and with no need for 
charge self-consistency furthermore. 

(2) Using the obtained wave functions, we can calculate the inner products among them 

and evaluate the N x N matrix = {m, k^^i, ky\n, k^^i^i, ky), where m and n are band 

indices of occupied states ranging from 1 to A^. The F matrix®^ is nothing but the discretized 
version of Berry connection A^. 

(3) Now we can construct the NxN matrix D as the product of F, D{ky) = FQ^iFi^2---FNk-2,Nk-iFNk-ifi- 

(4) Diagonalization of the D{ky) matrix gives us N eigenvalues \n{ky), whose phase angle 
can be obtained as 9n{ky)=lm[\ogXn{ky)]. 

This phase angle 9n{ky), which is now band dependent, is exactly what we need to 
determine the winding number and the corresponding topological numbers, as discussed in 
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FIG. 8: (Color online) Example of Wilson loop method for the 3D strong topological 
insulator TIN with index (1;000). The calculated 6m{ky) for all occupied bands [m is 
the band index) are plotted as functions of ky. The dashed lines are arbitrary reference 
lines. The Om{ky) crosses the reference line an odd number of times in half of ky direction 
{ky G [0,7r] or [tt, 27r]) for (a) k^ = 0, (c) ky — 0 and (e) kz = 0 plane (left panel), but an 
even number of times for (b) k^ = tt, (d) ky = tt and (f) kz = tt plane (right panel). 

Figures from Ref. 128. 


section II. In practice, for a system with many bands, we can draw an arbitrary reference line 
on the cylinder surface, parallel to the ky axis, and then count how many times the evolution 
lines of dn{ky) cross the reference line®^. For the calculation of Z 2 in time-reversal invariant 
systems, we can further reduce the computational task by considering only the Q < ky < tt 
part. As an example, we show in Fig. 8 the calculated 9n and their evolutions for one 3D 
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topological insulators TIN proposed recently^^®. All the six planes, each of them passing 
through four TRIM points, are shown. We see that the 9n lines cross the reference line an 
odd number of times for all three planes containing F, but an even number of times for those 
planes containing (7r,7r,7r). Therefore, we conclude that it is topologically non-trivial with 
Z 2 index (1;000). 

Finally, through the whole calculation process, we can always take advantage of Wannier 
functions, particularly the maximally localized Wannier functions (MLWF) as proposed by 
Vanderbilt^^®’^^®. For example, after the self-consistent electronic structure calculations, we 
can construct the MLWF for the energy bands close to the Fermi level, which can usually 
be reproduced up to a very high accuracy. In this process, however, some of the unoccupied 
bands should be included^^®. Once the MLWF are obtained, we can easily calculate the 
eigen wave functions and their inner products for any k-points using the interpolation tech¬ 
nique^^’ This will greatly simplify our task for the computation of topological numbers. 
In addition, the MLWF are also useful for the calculations of boundary states, as discussed 
below. 

III.C. Boundary States Calculations 

After the determination of bulk topological numbers by either of the methods presented 
above, the calculation of boundary states of real compounds becomes necessary. This is not 
just for the purpose of confirming topological properties - it is also very useful for experi¬ 
ments, such as angle-resolved photo-emission spectroscopy (ARPES) and scanning tunneling 
microscope (STM), because it provides necessary information for direct comparison. 

In practice, two methods are usually used for the calculation of boundary states. First, 
the boundary states can be directly calculated from first-principles calculations by using a 
supercell of slab geometry, as is usually done for other surface or edge calculations^^^. In 
this way, the atomic details, such as the surface or edge terminations, the atomic relaxations 
or reconstructions, the possible defects or absorptions, can be carefully treated and studied. 
This method is straightforward, although it is computationally demanding. Due to the 
presence of boundary states, we have to use a very large supercell with a thick vacuum 
region to avoid possible mutual coupling between the boundaries. The convergence test 
as a function of cell size is usually necessary. This method can generate detailed surface 
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eigen states, on the basis of which the spin texture of topological surface states (showing 
characteristic spin-momentum locking effect) can be obtained. With spin-orbit coupling 



included, the eigen states are expressed as two-component spinors, ipnk = I I , anc 

y 0nk4, J 

the expectation value of the spin operator S = |cr can be calculated easily by 


In addition, the layer-dependent distribution of surface states (i.e, the penetration depth) 
can also be studied carefully through these eigen states^°^. 

Second, we can take advantage of MLWF and calculate the boundary states from the 


Green’s functions of a semi-infinite system^^^. We can use the MLWF as a basis to con¬ 


struct an effective low-energy Hamiltonian Hwf, which can be regarded as a tight-binding 
model with its hopping parameters determined from the first-principles electronic structure 
calculations. The maximally localized property of MLWF guarantees the short range hop¬ 
ping. Using this Hwf as building blocks, we can construct a large Hamiltonian Hgemi for 
the semi-infinite system with only one boundary. This Hgemi is a block tridiagonal matrix 
if we consider only the hopping terms between the nearest-neighbor building blocks. In 
practice, we find that the thickness of the building block should also be tested to ensure 
convergence since in some cases the hopping among MLWFs might extend to more than 
the nearest-neighbor. Then we can use the iterative method^^^^^^® to solve the problem 
and get the projected Green’s functions onto the boundary. From these boundary Green’s 
functions, we can get necessary information such as the charge density of states and the spin 
density of states at the boundary^^®’^^^. This procedure is an approximate way and it cannot 
treat the atomic details at the boundary precisely, however, it provides the most important 
information for the gapless nature of the boundary states of topological insulators. 

IV. MATERIAL PREDICTIONS AND REALIZATIONS OF QAHE 

The topological electronic states discussed in previous sections are important and inter¬ 
esting, not only because they are conceptually new but also because they are realizable and 
testable in condensed materials. In recent years, this field has been developing very fast. 
Many candidate compounds with non-trivial topological electronic states are proposed, and 
some of them have been successfully synthesized and confirmed experimentally. The re¬ 
alizations of various topological electronic states have strongly promoted this field, and a 
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lot of interesting experiments on their non-trivial topological properties are now becoming 
possible. In this section we will focus on the topic of how to realize the QAHE in real 
materials. 

In spite of great success in recent studies of TIs and QSHE, the realization of QAHE 
has taken a long time. The key to achieving this goal is to find a proper Chern insulator. 
Any realistic material system in which QAHE can be realized must satisfy the following 
conditions: 

1. It must be two dimensional; 

2. It must be insulating in the bulk; 

3. It must break the TRS with a certain magnetic ordering (the simplest case is a EM 
long-range order); 

4. The occupied bands must carry a non-zero Chern number (therefore strong SOC is 
usually required). 

It is not difficult for a material to satisfy one or two conditions in the above list. But 
finding a realistic material that satisfies all of these conditions simultaneously turns out 
to be very challenging. This is the reason the Chern insulator is theoretically proposed 
long before the discovery of TIs, but its realization came much later. Nevertheless, the 
discovery of QSHE and have greatly stimulated the field of QAHE. It is 

natural to expect that starting from known 2D TIs and then breaking the TRS will be one 
of the simple ways to achieve the QAHE. The 2D TIs already satisfy three out of those four 
conditions, and these TIs can be effectively viewed as two Chern insulator layers related 
by TRS. Breaking of TRS can, in principle, destroy one layer and keep the other layer 
active, which will lead to the QAHE. Based on this idea, we can go further: if the exchange- 
splitting is strong enough, we can even realize the QAHE starting from a 2D narrow gap 
semiconductor, where the strong exchange-splitting may cause band inversion for half of the 
states (i.e, for only one of the spin subspaces, as will be addressed below). Unfortunately, it 
is still very hard to induce exchange-splitting in 2D TIs or narrow gap semiconductors; two 
known procedures, the magnetic proximity effect and magnetic doping, are most commonly 
employed. Despite the issue of magnetic ordering temperature Tc (which can be very low), 
magnetic doping will obviously introduce impurities that may destroy the insulating state 
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of parent compounds. On the other hand, using the magnetic proximity effect may have the 
advantage of avoiding the problem of impurity doping, while the possible exchange-splitting 
induced by this procedure is typically very weak. Therefore, in practice, we have to choose 
the candidate compounds and procedures carefully in order to balance various conditions 
and requirements. In this sense, the first-principles calculations play important, predictive 
roles for the final realization of the QAHE. As listed below, up to now, there have been 
several different types of proposals for the realization of the QAHE, and most of them start 
from TIs as the parent compounds. 

• To sandwich a TI thin film with EM insulators. The EM proximity effect will open 
up a gap for the Dirac surface states on both surfaces of the TI, which gives rise to 
the QAHE^h 

• To dope the magnetic ions into the TI thin film and make it a EM semiconductor 

The QAHE can be obtained by proper control of the doping. This is the only one of 
these proposals that has succeeded so far^^^’^^^. 

• To make a thin film (or quantum well structure) of a EM Weyl semimetal (such 
as HgCr 2 Se 4 ). The quantum size effect gives rise to the quantization of the crystal 
momentum along the growing axis, which at the same time quantizes the anomalous 
Hall coefficient at a certain film thickness^®’^^^. 

• To grow graphene, silicene, or other honeycomb lattice on top of magnetic insulators 
or with magnetic adsorbates. The low energy bands (which are tunable by external 
electric field) acquire both the Zeeman exchange-splitting and the SOC, which drive 
the system towards QAHE^^^^^®^. 

• To mimic the honeycomb lattice by growing double-layer thin films of a transition 
metal oxide with perovskite structure along the (111) direction, or by forming a sheet 
of organic moleculars that contain transition-metal ions. The magnetic long-range 
order might be stabilized among the d electrons of the transition metal. The QAHE 
can be realized in such systems if a suitable exchange field is established^^®"^®^. 

• Other proposals, such as heavy metal on magnetic insulator substrate^®^, strained 
epitaxial film with interface band inversion in EuO/GdN^®^ or CdO/EuO^®®, etc. 


42 


m<0 nn>0, without SOC m>0, with SOC 

FIG. 9: (Color online) Illustration of the band inversion mechanism. Schematic plot for 
the band structure with normal order (m<0), inverted order (m>0) without SOC, and 
inverted order with SOC, respectively. (See main text for detailed explanations). 

Most of the proposals mentioned above are supported by first-principles calculations, 
and we will start this section with a brief introduction about how to explore new topological 
materials from the viewpoint of a band inversion mechanism. Then we will discuss details of 
some of the proposals for the prediction and realization of the QAHE, with no attempting 
to cover them all. 

IV.A. Band Inversion Mechanism 

Though we have found, in the previous sections, that the Wilson loop method can reduce 
computation effort in calculating various topological invariants, it is still quite tedious and 
inefficient for wide-range searching or design of topological materials. From this point of 
view, band inversion is the most intuitive picture and most practically useful guideline for 
the initial screening of possible topological materials. Band inversion involves the energy 
order switching of low energy electronic bands around a certain TRIM point in the Brillouin 
Zone when compared with the energy ordering in the atomic limit (which is topologically 
trivial). Band inversion may happen when atoms form crystalline solids and acquire strong 
enough band dispersions. The mechanism was first demonstrated for a HgTe/CdTe quantum 
well structure^^ and later was found to be applicable to many other materials, including 2D 
and 3D Chern insulators^^®, topological Weyl semimetals^^’^® and topological 

Dirac semimetals^®’®^. 
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As schematically illustrated in Fig. 9, we consider two bands close to the Fermi level 
around the F point. Far away from the F point, the red state is occupied and energetically 
lower than the blue state. However, close to F point, they have an inverted energy ordering 
with the blue state lower and occupied at F. Let us assume that two bands do not couple 
when the SOC is absent, and they must cross at certain fc-point away from the F point. For 
time-reversal-invariant systems with SOC, the spin degree of freedom has to be taken into 
account, and at least four bands must be included in this simple model. The band crossing 
points are, in general, not stable and will open up a gap (in the presence of SOC). As a 
result, the system becomes an insulator if the chemical potential is located within the gap, 
and most importantly, this gapped state may be topologically non-trivial depending on the 
coupling terms. It should be noted that if additional crystal symmetries are considered, the 
four-fold degenerated band crossings are possibly protected. This situation is demonstrated 
in topological Dirac semimetals®®’®^. 

To be more explicit, we consider two parabolic bands with energy dispersion F'±(k) = 
±(—m -|- k^). The -|- and — signs are for the blue and red bands, respectively. The band 
inversion (around k=0) happens for the parameter region of m > 0. We first focus on a 2D 
system with broken TRS. Using two bands as a basis, we can construct an effective 2x2 
Hamiltonian, 


^e//(k) 


E+{k) M{k)* 

M(k) U_(k) 


(42) 


where M(k) is the coupling between the two bands that determines the topological properties 
of the resulting electronic structure. If m < 0, i.e., without band inversion, we will in general 
expect nothing but a normal insulator. However, if m > 0 for the band-inverted system, 
we may find many interesting topological phases, depending on the choice of M(k). For 
example, if M(k) = (k^ ± iky)^ {n > 1), we will get a Chern insulator with Chern number 
Z = n (for a 2D system). For a 2D system with TRS, the two bands should be replaced 
by four bands with each color (blue and red) representing a Kramers pair of bands, and 
the model Hamiltonian becomes 4x4 (block diagonal in its simplest form, such as the BHZ 
Hamiltonian^^). If we now have a form of M(k) = {k^ ± ikyY'^~^, we will expect a 2D 
topological insulator with ^2=1. In the presence of SOC, these forms of coupling term 
M(k) are in general allowed, and they depend on the material’s details. Similar discussions 
can be extended to 3D systems, and in this case, we can expect either 3D TIs^^° or 3D 
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topological semimetals^®’^®’®^. In real compounds, the band inversion may happen between 
s-p, p-p, p-d, or d-f states, which leads to different material classes of TIs^®^. 

The band inversion mechanism discussed above is practically useful, yet a crude picture 
for real compounds. First, this picture comes from the low energy effective Hamiltonian of a 
continuous model. If we put the model onto a lattice, we have to be careful about the number 
of times band inversion occurs. Particularly for the Z 2 TIs, if band inversion happens an 
even number of times or it happens around an even number of TRIM points, we will return 
to a trivial insulator phase (in the sense of .^ 2 = 0 , although other topological phases, such as 
topological crystalline insulator^^^’^^^, may be defined). Second, real compounds may have 
many low energy bands that are strongly fixed and must be carefully analyzed. Third, the 
possible M (k) terms can also be complicated and need to be treated carefully. In any case, 
a final conclusion must be carefully made based first-principles calculations as presented in 
previous sections. 

IV.B. QAHE in Magnetic Topological Insulators 

Several realistic systems, including Mn doped HgTe quantum wells^^®, magnetic-impurity 
doped Bi 2 Se 3 thin films^^®, Mn doped InAs/GaSb quantum wells^^°, magnetically doped 
junction Quantum Wells^^^ etc, have been proposed to support the QAHE. The essence of 
these proposals is to break the TRS in a 2D slab of QSHI by magnetic doping. If the FM 
order can be established, the introduced exchange splitting may destroy the band inversion 
in half of the states (say, for one spin channel), while keeping the band inversion in the other 
half of the states, leading to the QAHE. 

A HgTe quantum well with an appropriate thickness realized between two CdTe barriers 
can achieve the 2D QSH state^^’^^. Based on this structure, Liu et ah proposed that the 
QAHE can be obtained by magnetic Mn-doping in HgTe/CdTe quantum wells^^*, if the 
doped Mn ions are ferromagnetically ordered with magnetic momentum polarized along 
the out-of-plane direction. Unfortunately, (Hgi_a;Mna;)Te is experimentally paramagnetic 
(rather than FM), which prohibits the realization of QAHE. Due to the paramagnetic nature, 
in practical measurement, a magnetic field is required to polarize the Mn magnetic moment. 
Therefore, it is hard to distinguish the Hall and the anomalous Hall contributions. With 
the same idea, Zhang et al. proposed that a QSHI with bulk band gap around 0.1 eV 
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can be realized in junction quantum wells comprising II-VI, III-V, or IV semiconductors, 
and proper magnetic ordering can drive the system into the QAH state. Experimentally 
again, we have to wait for the synthesis and the well-controlled doping of a sample. In the 
following, we will take Bi 2 Se 3 and Bi 2 Te 3 (well-known topological insulators), as paradigms 
to show how the QAHE state is achieved in their thin him form with Cr-doping. This is the 
only experimentally successful example so far^^^. 


(a) 




• A site 
A B site 
TCsite 





EIG. 10: (Color online) Crystal structure of the Bi 2 Se 3 family of compounds, (a) The 
hexagonal supercell containing 15 atomic layers, (b) Top view of a QL in the triangular 
lattice. Three sets of different sites, labeled as A, B and C sub-lattices, respectively, are 
presented, (c) The first BZ. Four non-equivalent TRIM points T (0, 0, 0), L(7r,0,0), 
F(7r,7r,0) and Z’(7r,7r,7r) are denoted in the 3D BZ. The corresponding surface 2D BZ is 
represented by the dashed hexagon, and f, M and K are the corresponding TRIM special 
k points in the surface BZ. Figure from Ref. 104. 

The Bi 2 Se 3 family of compounds have a rhombohedral crystal structure with space group 
(R3m). The system has a layered structure with five atomic layers as a basic unit, named 
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a quintuple layer (QL), as shown in Fig. 10. The materials Bi 2 Se 3 , Bi 2 Te 3 and Sb 2 Te 3 have 
been theoretically predicted and experimentally observed to be topological insulators with 
a large bulk band gap^^^®^. Recent experimental progress has shown that well-controlled 
layer-by-layer MBE thin film growth can be achieved^®and various magnetic transition 
metal elements (such as Ti, V, Cr, Fe) can be substituted into the Bi 2 Se 3 family of parent 
compounds with observable ferromagnetism even above 

We first discuss how to achieve a FM insulating phase in a semiconductor doped with 
dilute magnetic ions. To simplify the discussion, we divide the whole system into two sub¬ 
systems describing the local moments and the other band electrons, respectively. We further 
assume that the magnetic exchange among local moments is mediated by the band electrons. 
If we consider only the spatially homogeneous phase, the total free energy of the system in 
an external magnetic field H can be written as 

i - Jet/AhM, - (Afi + AQ H, (43) 


where XL/e is the spin susceptibility of the local moments/electrons, M^/e denotes the mag¬ 
netization for the local moment and electron sub-system and Jg// is the magnetic exchange 
coupling between them. The spontaneous FM phase can be realized when the minimization 
procedure of the free energy gives a non-zero magnetization in the zero magnetic field limit 
(i/=0). This leads to the requirement of Jeff~XL^X7^ > 0) equivalently Xe > ^/{JlffXL)- 
This result means that in order to have a non-zero FM transition temperature (Tg), a sizable 
Xg is needed. Within linear response theory, the spin susceptibility Xe for a insulator is given 
by the van Vleck formula 


„2 (^^1 Sz \ck) {ck\S, \vk) 


v,c,k 


^ck ^'1 


vk 


(44) 


where po is the vacuum permeability, pe is the Bohr magneton, \ck) and \vk) are the 
Bloch functions in the conduction and valence bands respectively, Eck and E^k are the eigen 
energies of \ck) and \vk), and Sz is the spin operator of electrons and we consider only the 
^(-direction, i.e. perpendicular to the 2D plane, for simplicity. 

In most dilute magnetic semiconductors, such as (Gai_a;Mn 2 ;)As, which has s-like conduc¬ 
tion and p-like valence bands, the matrix elements {ck \ Sz \vk) are very small. The electronic 
spin susceptibility is therefore negligible for the insulating phase, and a finite carrier density 
is required to mediate the magnetic coupling (for example, by the RKKY mechanism^^®). 
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FIG. 11: Van-Vleck type spin susceptibility of Bi 2 Se 3 bulk as a function of spin-orbit 
coupling (SOC) strength (Aq is the actual SOC strength). Figure from Ref. 139. 

For the Bi 2 Se 3 family of materials, however, the semiconductor gap is opened by the SOC 
between two p-orbital bands. The spin operator Sz has finite matrix elements between the 
valence and conduction bands, which are further enhanced due to the band inversion^^®. 

In Fig. 11, we show the calculated spin susceptibility of bulk Bi 2 Se 3 . At the F point, the 
band inversion occurs when the relative SOC strength A/Aq exceeds 0.5, and after that the 
spin susceptibility starts to increase rapidly. At the actual SOC strength, Bi 2 Se 3 -family 
materials have a considerable spin susceptibility xT- Hence the local moments of the mag¬ 
netic impurities in such materials can be ferromagnetically coupled through the van Vleck 
mechanism^^®. Recently the interplay between magnetism and the topological property in a 
material has been discussed in more details in Ref. 180 and 181. 

From first-principles calculations, we find that the insulating magnetic ground state dis¬ 
cussed above can indeed be obtained by a proper choice of magnetic dopants. Experiments^^^”^^^ 
and theoretical calculations^^®’^®^ show that the magnetic dopants will mostly substitute for 
Bi ions, which has a nominal 3-|- valence state. In order to avoid introducing free carriers 
into the parent material when we dope the system, it is natural to choose those transition 
metal elements that also have a stable 3-|- chemical state, such as Ti, V, Cr and Fe. Self- 
consistent first-principles calculations for Bi 2 Se 3 doped with Ti, V, Cr and Fe have been 
performed. The calculated densities of states (DOS) shown in Fig. 12 suggest that an 
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FIG. 12: Calculated densities of states (DOS) for Bi 2 Se 3 doped with different transition 
metal elements. The Fermi level is located at energy zero, and the positive and negative 
values of DOS are used for up and down spin, respectively. The lines are projected partial 
DOS of the 3d states of transition metal ions. It is shown that Or or Fe doping will give 
rise to the insulating magnetic state. Figure from Ref. 139. 


insulating magnetic state is obtained only for Or or Fe doping, while the states are metallic 
for Ti or V doping cases. 

The above results can be understood as illustrated in Fig. 13. We first point out that, 
for all the cases, the dopants are nearly in the 3+ valence state, and we always obtain 
the high-spin state due to the large Hund’s coupling of 3d transition metal ions. For the 
Fe-doped case, the Fe^"*" favors the d^^d^^ configuration in a high-spin state, resulting in a 
gap between the majority and minority spins, which directly leads to the insulating state 
in Fe-doped samples as shown in Fig. 13a. For the Cr-doped case, Cr ions substitute in 
the Bi sites and feel an octahedral crystal field formed by the six nearest neighboring Se^“ 
ions. Such a local crystal field splits the d-shell into t 2 g and Sg manifolds. This splitting 
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is large enough to stabilize the configuration of a Cr^'*' ion, resulting in a gap 

between the t 2 g and Eg manifolds, as shown in Fig. 13(b). For the case of Ti or V doping, 
the t 2 g manifold is partially occupied, leading to the metallic state, as shown in Fig. 13(c). 
We note that, although the local density approximation in density functional theory may 
underestimate the electron correlation effects, the inclusion of electron-electron interaction 
U (such as in the LDA-|-?7 method) should further enhance the gap (it may also reduce the 
p-d hybridization and Tc^®^). 

The exchange field A can be estimated as A = xJeff{S) at the mean field level, where 
X is the doping concentration and (S) is the mean field expectation value of the local spin, 
Jeff is the effective exchange coupling between the local moments and the band electrons, 
which is estimated to be around 2.7 eV for Cr-doping and 2.8 eV for Fe-doping in Bi 2 Se 3 
by first-principles calculations. With the concentration of the magnetic dopants at 10%, the 
FM Curie temperature can be estimated in the mean field approximation^®^ to reach the 
order of tens of K. 

Once the FM order is achieved in TI at a reasonable temperature, the QAHE can in 
principle be realized in 2D thin films of such systems by tuning the exchange splitting (i.e, 
the density of magnetic doping). Since the bulk states are always gapped, we focus on the 
simplest low energy effective Hamiltonian consisting of Dirac-type surface states only: 


H = H,f + H, 


sf “T -n. Zeeman 


—VF{k X a)z ml 

+ 

gMa^ 

0 

rufc vpik X a)z 

0 

gMa^ 


(45) 

with the basis of |t t); 1^ i) t) > 1^ 4'); where “t” (“6”) represents the surface states on the 

(a) (b) (c) 
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FIG. 13: The configuration of the d electrons for the (a) Fe, (b) Cr, (c) Ti and V doping 

cases. 


50 












(a) 




FIG. 14: Evolution of the subbands structure upon increasing the exchange field. The red 
color denotes the subbands that have even parity at the F point, and blue color denotes 
subbands with odd parity at the F point. The dashed lines denote the spin up electrons; 

solid lines, spin down electrons, (a) The initial subbands are not inverted. When the 
exchange field is strong enough, a pair of inverted subbands appears (red solid line and 
blue dashed line), (b) The initial subbands are already inverted. The exchange field 
releases the band inversion in one pair of subbands (red dashed line and blue solid line) 
and increases the band inversion in the other pair (red solid line and blue dashed line). 


top (bottom) surface, and “t”, “4-” represent the spin up and down states, respectively, vp 
is the Fermi velocity, ruk describes the tunneling effect between the top and bottom surface 
states, g is the effective g-factor, a are the Pauli matrices, and M represents the exchange 
field along the z-direction, introduced by the FM ordering. For simplicity, spatial inversion 
symmetry is assumed, which requires that vp, g and M take the same values for top and 
bottom surfaces. In thick enough slab geometry [rrik ~ 0), the spatially separated pairs of 
surface states are well defined for the top and bottom surfaces. However, with the reduction 
of the film thickness, quantum tunneling between the top and bottom surfaces becomes 
more and more pronounced, giving rise to a finite mass term which can be expanded 
as nik = rriQ + B{k‘l + ky) up to the second order of k. The Hamiltonian can be rewritten 
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in terms of the symmetric and anti-symmetric combination of the surface states on top and 
bottom surfaces as 


H = Hsf + H. 


^Zeeman 


hk + gMa^ 0 

0 hi- gMa^ 


(46) 


with the following new basis: |-|- t), |— 4-)) 1+ 4-)) |~ t); where |± t) = (|^ t) i 1^ t))/'s/2, 
|± 4 -) = {\t 4 -) ± 1 ^ 4 -)) /V^- Here h{k) = mkCTz + vp {kyGx — k^ay), which is similar to the 
BHZ model describing the low energy physics in HgTe/CdTe quantum wells^^. When rrioB < 
0, band inversion occurs, and the system will be in the QSH phase if this is the only band 
inversion between two subbands with opposite parity. Regardless of whether this condition 
is satisfied or not, a strong enough exchange field will induce the QAH effect in this system, 
thanks to the presence of the matrix in the exchange field [gMa^) and the opposite signs 
of the Zeeman coupling terms in the upper and lower blocks of the effective Hamiltonian 
in Eq. (46). The exchange field breaks the TRS and increases the mass term of the upper 
block, while reducing the mass term for the lower block. More importantly, a sufficiently 
large exchange field can change the Chern number of one of the two blocks. There are 
two cases, as illustrated in Fig. 14. In the first case, the four-band system is originally in 
the topologically trivial phase; the exchange field will induce band inversion in the lower 
block and push the two sub-bands in the upper block even farther away from each other. 
Therefore, the 2D model with a negative mass in the lower block contributes —e^/h for 
the Hall conductance (see Fig. 14(a)). In the other case, the system is originally in the 
topologically non-trivial phase, and both blocks have inverted band structures. A sufficiently 
large exchange field can increase the band inversion in the lower block and release the band 
inversion in the upper block. Again the negative mass in the lower block contributes —e ^/h 
for the Hall conductance (see Fig. 14(b)). Such a mechanism is general for the thin films of 
TIs with FM ordering, it is guaranteed by the fact that the surface states on the top and 
bottom surfaces have the same p-factor. 

In Fig. 15 A-C, we plot the lowest four sub-band levels at the T point as functions of 
exchange field. Level crossings between the lowest conduction bands (blue lines) and valence 
bands (red lines) are found for all three values of layer thickness, signaling the quantum phase 
transition to the QAH state. In the insulator case, where the chemical potential is located 
inside the energy gap between the conduction and valence sub-bands, the Hall conductance 
is determined by the Chern number of the occupied bands and must be an exact integer 
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FIG. 15: Quantized anomalous Hall (QAH) conductance. (A), (B) and (C) are the lowest 
sub-bands at the F point, plotted versus the exchange held, for Bi 2 Se 3 hlms with 
thicknesses of 3, 4 and 5 quintuple layers (QL), respectively. The red lines denote the top 
occupied (bottom unoccupied) states. With an increase of the exchange held, a level 
crossing occurs (indicated by arrows), indicating a quantum phase transition to QAH 
state. (D) Calculated Hall conductance for 3, 4 and 5 QL Bi 2 Se 3 hlms under ferromagnetic 
exchange held. As we expect, the Hall conductance is zero before the QAH transition, but 

Wy = e^//i afterward. Figure from Ref. 139. 


multiple of e^//i. The calculated Hall conductance for Bi 2 Se 3 thin hlms of three typical 
thicknesses (Fig. 15 D) jumps from 0 to 1 at the corresponding critical exchange held for 
the level crossing. The QAH effect with higher Hall conductance plateaus in the Cr-doped 
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FIG. 16: The sub-bands dispersion of 2D Bi 2 Ses thin film are plotted for SQL slab with 
exchange field (A) -0.05 eV, (C) -0.24 eV. The corresponding Hall conductance are plotted 
in the right panels as the function of the chemical potential. The zero chemical potential 
corresponds to the ideal stoichiometry case. The systems in (A) is normal insulator, which 
has zero platforms in Hall conductance when the chemical potential locates in the energy 
gap. The system in (C) is QAH insulators, in which the Hall conductance has quantized 

value e^//i . Figure from Ref. 139. 


Bi 2 (Se,Te )3 was further discussed in Ref. 185, where the stronger exchange field leads to 
more band crossings and a larger Chern number. 

In Fig. 16, we plot the sub-band dispersion and corresponding Hall conductance as 
functions of chemical potential for 3 QL Bi 2 Se 3 thin films with different exchange fields. A 
quantized plateau in Hall conductance should be observed when the chemical potential is 
inside the gap. For the QAH phases (Fig. 16 (D)), a non-zero integer plateau is observed 
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when the chemical potential falls into the energy gap; while for the non-QAH phases (Fig. 
16 (B)), the corresponding Hall conductance is zero. In a recent experiment, Chang, et al.^^^ 
prepared the Cr-doped (Bi,Sb) 2 Te 3 thin films with well controlled chemical potential and 
long-range ferromagnetic order. Careful control of the Bi and Sb concentration is needed to 
separate the surface states from the bulk states. By tuning the gate voltage (and therefore, 
the chemical potential), they observed the quantization of the Hall resistance at h/e^ at zero 
external magnetic filed with temperature around 30 mK, accompanied by a considerable drop 
in the longitudinal resistance. A similar experiment was further performed by Checkelsky, et 
al.^^^, confirming the existence of QAHE in this system. For readers who want to know more 
experimental details, please refer to their original papers. To conclude this part, we finally 
point out that the QAHE can be induced even by in-plane magnetization, as discussed in 
Ref. 186. 

IV.C. QAHE in thin film of Weyl semimetals 

As discussed in sec. H.C, Weyl semimetal is a new topological state of 3D quantum 
matter, unlike the 3D topological insulators. It can be characterized by Weyl nodes at the 
Fermi level in the bulk and Fermi arcs on surfaces. Around Weyl nodes, the low-energy 
physical behavior is like that of 3D two-component Weyl fermions®® and can be described 
by the Hamiltonian H = Vfcr where Vf is the Fermi velocity, cr is the vector of the Pauli 
matrices, and k is the momentum as measured from the Weyl node. A Weyl fermion is 
half of a Dirac fermion, and Weyl fermions must appear in pairs. The Weyl Hamiltonian is 
robust against perturbations since it uses all three of the Pauli matrices. The perturbations 
only shift the position of the gap-closing point in momentum space and the Weyl nodes 
can be removed only when a pair of them meet in momentum space. To obtain a Weyl 
semimetal, either time-reversal symmetry or inversion symmetry must be broken, and there 
have been several proposals to realize In particular, Weyl semimetals are 

closely related to the chiral anomalyand the QAHE can naturally be expected in the 
quantum well structure of Weyl semimetals with broken TRS. In this section, based on first- 
principles calculations, we will show that HgCr 2 Se 4 is a EM Weyl semimetal with Eermi arcs 
on its surface, and its 2D thin film is a promising system for realizing the QAHE. 

HgCr 2 Se 4 is a FM spinel exhibiting large coupling effects between electronic and magnetic 
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FIG. 17: (Color online) Crystal structure of HgCr 2 Se 4 spinels. 


properties^®®. The spinel structure (space group Fd3m) of HgCr 2 Se 4 with Cr®+ octahedrally 
surrounded by Se ions, as shown in Fig. 17, yields a half-filled t 2 g shell and a spin number of 
S=3/2. First-principles calculations^® confirm that FM solution is considerably (2.8 eV/f.u.) 
more stable than non-magnetic solution, and the calculated moment (6.0 /i^/hu.) is in good 
agreement with experiments^®®’^®®. Without SOC, the electronic structures shown in Fig. 
18(a) and (b) suggest that it is nearly a “zero-gap half-metal”. It is almost a half-metal 
because a gap exists only in the up-spin channel just above the Fermi level; it is nearly 
zero-gapped because of the band-touching around the F point just below Fermi level in 
the down-spin channel. Quite large exchange splitting and strong octahedral crystal field 
together lead to a high spin state of Cr®"*" ions in a configuration with a gap 

between the and manifolds as schematically shown in Fig. 18(c). The Se-4p states 
(located from about -6 eV to 0 eV) are almost fully occupied and contribute to the top part of 
the valence band dominantly. The hybridization with Cr-Sd states makes the Se-4p slightly 
spin-polarized with the opposite moment (about -0.08 /r^/Se). The zero-gap behavior in 
the down spin channel is the most important characteristic here (Fig. 18(b)), because it 
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suggests the inverted band structure around F, similar to the case in HgSe or HgTe 
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FIG. 18: (Color online) Electronic structures of HgCr 2 Se 4 . (a) The total and projected 
partial density of states (DOS); (b) The band structures without SOC (showing the up and 
down spin parts separately); (c) Schematic plot of band-inversion, where the |S', j.) state is 
lower than the \P) states; (d) The band structure including SOC (with majority spin along 
the (001) direction). The detail around F is enlarged in inset. Figure from Ref. 56. 


Taking the four (eight if spin is considered) low energy states at the F point as basis, 
which are identihed as \Px), \Py), \Pz), and |S'), the low energy physics becomes the same 
as in HgSe or HgTe, and the only difference is the presence of exchange splitting. The band 
inversion (see Fig. 18(c); [S', j,) is lower than |P)) is due to the following two factors: 1) 
Hg-5d states around -7.0 eV below Fermi level are very shallow and their hybridization with 
Se-4p states has pushed the anti-bonding Se-4p states quite high, similar to the situation 
in HgSe. 2) The hybridization between unoccupied Cr-Sd^ and Hg-6s'*' states will push the 
Hg-6s'*' state to lower in energy. These two key factors lead to band inversion, with the l^, 4-) 
lower than the \P3, |) states by about 0.4 eV. This value is further enhanced to be 0.55 eV 
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FIG. 19: (Color online) Weyl nodes and gauge flux in HgCr 2 Se 4 . (a) The band-crossing 
points in fc-space; (b) Chern number as function of kz\ (c) Schematic plot of the band 
dispersion around the Weyl nodes in the kz=:tkl plane, and the inset shows the chiral spin 
texture, (d) Gauge flux evaluated as Berry curvature in the [k^^ k^) plane. Figure from 

Ref. 56. 


in the presence of SOC. Further LDA+G calculations (with effective U around S.OeV^®^’^®^) 
indicate that the band-inversion remains even if the correlation effect is properly considered, 
unless U is unreasonably large (> 8.0 eV). Experimental observations of metallic behavior 
at low temperature for various samples^®^”^®^ strongly support this conclusion of inverted 
band structure. 

In the presence of SOC, the low energy eigenstates at F are given as ||,±|), ||,±i), 
||,±i), and IS”, ±|), which can be constructed from the |P) and \S) states, similar to the 
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case of HgSe and HgTe again. Considering the exchange splitting, the eight states at F are 
separated in energy. The ||, |) has the highest energy, and the — ^) has the lowest. Due 
to the band-inversion, several band-crossings are observed in the band structure. Among 
them, however, only two kinds of band-crossings (called A and B) are important for the 
states very close to the Fermi level, as shown in Fig. 18(d). As schematically shown in Fig. 
19(a), the A-type crossing gives two points located at = ±A:^ along the V — Z line. For 
the 2D planes with fixed-A:^ [k^ ^ 0, ±A:^), the in-plane band structures are all gapped in the 
sense that a curved Fermi level is defined. Therefore, the Chern number C for each fc^-fixed 
plane can be evaluated. It is found that the crossing A (i.e., \kz\ = kf) is located at the phase 
boundary between the C = 2 and C = 0 planes, i.e., (7 = 0 for the planes with \kz\ > k^, 
while (7 = 2 for the planes with \kz\ < k'^ and kz ^ 0- This change in Chern number 
indicates that these two crossing points (Weyl points) are topologically protected so long 
as the Weyl points remain separated in momentum space. The B-type crossing is a closed 
node line surrounding the F point in the kz=0 plane. It is just accidental and is due to the 
presence of crystal mirror symmetry with respect to the kz=0 plane. The B-type crossing 
is not as stable as the A-type crossing, in the sense that it can be avoided by changing of 
the crystal symmetry (such as the reorientation of magnetic momentum away from the z 
direction). 

Assuming that the magnetic momentum of the system is oriented to the z axis, the low 
energy physics here can be caught by the following effective two-band modeF®’^®®: 


Heff = 


M Dkzkl 
Dkzkl -M 


(47) 


in bases ||, |) and \S, —\)- Here k± = kx iky, and M = Mq — /5/c^ is the mass term 
expanded up to the second order of k, with parameters Mq > 0 and /3 > 0 to ensure band 
inversion. Since the two bases have opposite parity, the off-diagonal element has to be an 
odd function of k. In addition, k"^ has to appear to conserve the angular moment along z- 
direction. Therefore, to the leading order, kzk'^ is the only possible form for the off-diagonal 
element. The energy dispersion E{k) = -|- D‘^kl{k‘^ + kyY has two gapless solutions: 

one is the degenerate points along T — Z with kz = Ek^ = ±y/Mo//5; the other is a circle 
around the F point in the = 0 plane determined by the equation k^ + ky = Mq/(3. They 
are exactly the band-crossings obtained from the first-principles calculations. Due to the 
presence of k"^ in the off-diagonal element^®, it is easy to find that Chern number C equals 
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FIG. 20: (Color online) Edge states and Fermi arcs of HgCr 2 Se 4 . (a) Calculated edge 
states for the plane with A:2=0.067r. A ribbon with two edges is used, and there are two 
edge states for each edge since C=2. The states located at different edges are indicated by 
different colors, (b) Calculated Fermi arcs for the side {ky,kz) surface. Figure from Ref. 56. 


to 2 for the planes with \kz\ < k^ and kz ^ 0, otherwise C=0. The band dispersions near 
the Weyl nodes at kz = plane are thus quadratic rather than linear, with chiral in¬ 

plane spin texture, as shown in the inset of Fig. 19(c). The two Weyl nodes located at 
±kz have opposite chirality and form a single pair of magnetic monopoles. The gauge flux 
starting from and ending at them in A:-space is shown in Fig. 19(d). The band-crossing loop, 
i.e. node-line, in the kz = 0 plane is not topologically unavoidable; however, its existence 
requires that the gauge flux in the kz=0 plane (except the loop itself) must be zero. 

This Weyl semimetal (also called Chern semimetal) state realized in HgCr 2 Se 4 will lead 
to novel physical consequences, which can be measured experimentally. First, each fc^-flxed 
plane with non-zero Chern number can be regarded as a 2D Chern insulator, and there 
must be chiral edge states for such plane if an edge is created. The number of edge states 
is two for the case of C=2 (see Fig. 20(a)), or zero for the case of C=0. If the chemical 
potential is located within the gap, only the chiral edge states can contribute to the Fermi 
surface, which are isolated points for each plane with (7=2 but nonexistent for the plane 
with (7=0. Therefore, the trajectory of such points in the {k^, kz) surface or {ky, kz) surface 
form non-closed Fermi arcs, which can be measured by ARPES. As shown in Fig. 20(b), 
the Fermi arcs connect the Weyl points at kz=^kl, and are interrupted by the kz=Q plane. 
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This is very different from conventional metals, in which the Fermi surfaces must be either 
closed or interrupted by the Brillouin zone boundary. The possible Fermi arcs have also 
been discussed for pyrochlore iridates®^. 

The QAHE, on the other hand, is a unique physical consequence characteristic of the 

Chern semimetal nature of HgCr 2 Se 4 by considering its quantum-well structure. For 2D 

2 

Chern insulators, the transverse Hall conductance should be quantized as a^y = C^, where 
C is the Chern number. Considering the fc^-fixed planes, the Chern number C is non¬ 
zero for limited regions of k^, and this is due to the band inversion around F as discussed 
above. In the quantum well structure, however, those low energy states around F should be 
further quantized into subbands (labeled as |i7„) and \En) for hole and electron subbands, 
respectively), whose energy levels change as functions of him thickness. As shown in Fig. 
21(a), when the film is thin enough, the band inversion in the bulk band structure will be 
removed entirely by the finite size effect. With the increment of film thickness, finite size 
effect gets weaker and the band inversion among these subbands is restored subsequently, 
which leads to jumps in the Chern number or the Hall coefficient As shown in Fig. 

21(b), if the film is thinner than 2lA, the a^y is zero; once the film thickness is larger than the 
critical thickness, we find subsequent jumps of a^y in units of 2e^/h. In fact, the anomalous 
Hall effect has been observed for the bulk samples of HgCr 2 Se 4 ^®®, and more recently, n-type 
HgCr 2 Se 4 was confirmed to be a single s-band half metaF°°. This is in sharp contrast to 
pyrochlore iridates, in which the anomalous Hall effect should be vanishing due to the AF 
ordering^^. 

This proposal of realizing QAHE in thin films of FM Weyl semimetals has some advan¬ 
tages compared to the previous one for the magnetic-ion doped TIs thin films. First, it is 
based on the stoichiometric material and avoids the difficulty of doping, which is in general 
good for raising the mobility of the sample. Second, FM ordering has been established in 
bulk HgCr 2 Se 4 at temperature as high as 100 which makes QAHE plausible at high 

temperatures. Third, the proposed system can reach a QAHE with a high Chern number 
(i.e., Chern number 2 or even higher) and higher plateau^®^. Therefore, further experimental 
efforts in realizing QAHE in thin film of HgCr 2 Se 4 are much anticipated and very promising. 

IV.D. QAHE on Honeycomb Lattice 
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FIG. 21 : (Color online) Quantized Anomalous Hall effect in HgCr 2 Se 4 thin film, (a) The 
subband energy levels at F point as function of film thickness, (b) The Hall conductance as 

function of film thickness. Figure from Ref. 56. 


Since honeycomb lattice plays a unique role in exploiting non-trivial topology^^’^®, it is 
instructive to take a look at the geometric property of its electronic band structures, paying 
attention to possible configurations of Berry curvatures. There are two sites in a unit 
cell of honeycomb lattice, and with the nearest neighbor hoppings, the valance band and 
conduction band touch linearly and form the Dirac cones at the six corners of the Brillouin 
zone. Taking into account the C 3 symmetry, one can consider the two inequivalent Dirac 
cones at K, K' = (±1,0) with lattice constant d, which are called valleys and denoted by 

Tz = ±1. The Hamiltonian for low-energy physics is given by 

( 0 Tzkx + iky 

y 

Tzkx — iky 0 

where the momenta are measured from the valleys. This is the 2 D version of Eq. (35) with 
the basis of the 2 x 2 matrix referring to AB sublattices of honeycomb lattice. Therefore, the 
electrons are chiral at the valleys, which makes the honeycomb lattice special in generating 
non-trivial topology. 

Let us introduce a staggered field U on the AB sublattices, and rewrite the Hamiltonian 
(48) as 

Htz T~zkx(^X T kyCTy U<Tzi (49) 

with Pauli matrices cFx,y,z- The dispersion is given hy E = + k"^, with a gap opened 
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FIG. 22: (Colour online) Possible configurations of Berry curvatures on honeycomb lattice: 
(a) CDW, (b) QSHE, (c) SDW, (d) QAHE, (e) spin resolved QAHE characterized by 
simultaneous non-zero charge and spin Chern numbers^^h Red/blue for spin up/down 
channels. Edge currents are displayed in insets in (b), (d) and (e) for topologically 
non-trivial states. and are Chern numbers in spin-up and -down channels, and 
C'c/s = ± Cj, are charge and spin Chern numbers respectively. Eigure from Ref. 161. 

by the staggered potential U. Eor this insulating state, one can define a fictitious magnetic 
field d = {Tzkxiky^U) Hamiltonian (49) on the pseudo spin cr, which exhibits a topological 
structure around each valley known as meron. Its topological property can be character¬ 
ized by counting the winding number of the meron structure (i.e., the integration of Berry 
curvature over a patch covering the region of the meron, or in other words, the Berry flux 
passing through the 2D BZ in the vicinity of the each valley, divided by 27r) as^® 

^ = ^ j dkxdkyYi ■ (dfc^n x d^^n) = sgn(T^[/)/2, (50) 

with n = d/|d|. 

IV.D.l Combinations of Berry curvatures on honeycomb lattice 
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For a staggered electric potential i/eiec = y<^z where 1/ is a constant, the signs of the Berry 
curvature around two valleys are opposite, as seen in Eq. (50), which gives a topologically 
trivial insulator called a charge density wave (CDW) as displayed in Fig. 22(a). On the 
other hand, if a staggered magnetic flux is introduced to spinless electrons on honeycomb 
lattice^^, the staggered potential term can be written as (around two valleys), 

which depends on the valley variable in contrast to the staggered electric field. Since the 
topological charges are now given by c = sgn(t)/2, the Berry curvatures around two valleys 
are aligned, which characterizes a QAHE state (see Fig. 22(d) where spin is included in 
addition to the Haldane model). This staggered held can be generated by shining circularly 
polarized light on graphene, as revealed by the Floquet theorem^®^. Taking into account 
the true spin degree of freedom, the SOC term of i^so = can also be included as 

suggested in the Kane-Mele modeF®, and the SOC term again takes the form of a staggered 
field, and depends additionally on spin variable s^. In this case, the Berry curvatures around 
two valleys in the spin up and down channels are aligned individually and oppositely for 
different spin channels, resulting in the configuration in Fig. 22(b), which characterizes the 
QSHE state^®. Antiferromagnetic (AFM) exchange field i^afm = MszCTz is another possible 
staggered field, which reverses the sign of the Berry curvature in the spin-down channel for 
both valleys (compared to the CDW state shown in Fig. 22(a), assuming the same sign for 
M and V), resulting in the configuration in Fig. 22(c) associated with a topologically trivial 
state called spin density wave (SDW). Coupling of the valley degree of freedom to the AFM 
order in manganese chalcogenophosphates was recently discussed^°^. 

We observe that the configuration of Berry curvature given in Fig. 22(e) corresponds to 
an additional topological state^®^ (see also^°^), by flipping the sign of the Berry curvature 
around K' valley in the spin up channel from the CDW state. In contrast to QSHE with a 
non-zero spin Chern number Cg = — (7^^205,206^ which is equivalent to the Z 2 index, and 

QAHE with a non-zero charge Chern number Cc = this state is characterized by 

simultaneous non-zero charge and spin Chern numbers. We may call the state in Fig. 22(e) 
spin-resolved QAHE, as opposed to other QAHE states where spin information is lost. 

Note that configurations of the Berry curvature in Figs. 22(a), (b), (c) and (d) are con¬ 
verted to each other by flipping the sign of two Berry curvatures either at the same valley 
or with the same spin. This reflects certain symmetries in these states: TRS in Fig. 22(a), 
TRS and spacial IS in Fig. 22(b), IS in Fig. 22(d), and the combination of TRS and IS in 
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FIG. 23: (Color online) Band configurations for G = 0 (left) and |A — M\ < V < A + M 
(presuming A, M > 0) for (a) M < A and (b) M > A. Red (blue) curves are for spin-up 
(-down) channels, and A and B for the sublattices. The straight lines are for topological 
edge states in ribbon systems, and the dashed ellipses indicate the valleys where band 

inversions occur upon increasing V. 

Fig. 22(c) (although both symmetries are broken individually). All symmetries are broken in 
the state portrayed in Fig. 22(e). One can realize the state in Fig. 22(e) by using three helds 
simultaneously - SOC, AFM exchange field and the staggered electric potential - which 
makes fine control on electronic states possible: 

H = -\- kyGy -\- {XtzSz + Ms;, + V)a;,. (51) 

We reveal behaviors of the system upon tuning the electric potential V with the M field 
and A fixed. For A > M > 0, the system takes a QSHE state at R = 0 (see Fig. 22(b)). 
The degeneracy in band dispersions is protected by the product of TRS and IS, a pseudo 
TRS. The gaps are 2(A -|- M) and 2(A — M) at the two valleys, respectively. As illustrated 
in Fig. 23(a), when V increases from zero, the gap in the spin-up channel shrinks, closes at 
V = \ — M and reopens at the K' valley; this quantum phase transition makes the spin-up 
channel topologically trivial, while the spin-down channel remains topologically non-trivial. 
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This QAHE state is possible since the electric field V breaks the IS, therefore breaking the 
pseudo TRS. When V increases further, the gap in the spin-down channel shrinks, closes at 
V = A + M and reopens at K valley (not shown explicitly). This second quantum phase 
transition makes the spin-down channel trivial as well, and the system transforms into the 
topologically trivial CDW state. 

For M > A > 0, the system is in a topologically trivial SDW state at F = 0 (see 
Fig. 22(c)). As seen in Fig. 23(b), when V increases, the gap in the spin-down channel 
shrinks, closes aXV = M—X and reopens at the K' valley, which brings the spin-down channel 
into a topological state. When V increases further, the gap in the spin-down channel shrinks, 
closes at V = X + M and reopens at K valley, which drives the system to the topologically 
trivial CDW state. It is clear that the spin-resolved QAHF state is realized when the sizes 
(absolute values) of the three fields satisfy the condition of being able to form a triangle, 
which indicates clearly that the interplay among the three fields is crucial. In this case the 
sign of the total staggered field U = XtzSz + Msz -|- F at one valley in one spin channel 
is opposite to the other three cases, which yields the configuration of Berry curvature in 
Fig. 22(e). 

The above discussions for the possible topological phases in honeycomb lattice by fine 
tuning of parameters can be generally applied to many different material systems, such 
as graphene^^^^^^^, silicene^^^^^®^, perovskite bilayers^^®^^®^, etc. Valleys, as an additional 
degree of freedom of electrons on honeycomb lattice, have attracted considerable interest, 
for which valleytronics has been coined as a new concept^^^’^^® . We will leave the discussions 
of topological states in graphene and silicene to readers, and in the following, we concentrate 
on the perovskite bilayer systems, which are interesting and promising due to the presence 
of transition metal elements. 

IV.D.2 Honeycomb lattices by perovskite structures 

Material realization of topological insulators (TI) has been mainly focused on narrow 
band-gap semiconductors with heavy elements such as Hg or Bi, in which the electronic 
properties are dominated by s and p orbitals. Recently a completely different materials class 
- heterostructures of transition-metal oxides (TMOs) involving d electrons - was proposed^^®. 
The TMO heterostructures are readily available due to achievements in the fields of oxide 
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FIG. 24: (Color online) Perovskite structure as a buckled honeycomb lattice grown along 
the [ 111 ] crystal direction. Adapted from Ref. 156. 


superlattices and oxide electronics^°®^^^h The sample of layered structure of TMOs can be 
obtained with atomic precision^^^”^^^. This make the fine tuning of sample properties, such 
as the lattice constant, carrier density, SOC and electron-electron correlation, very plausible. 
Possible 2 D TIs have been proposed for bilayers of perovskite-type TMOs grown along the 
[ 111 ] direction of the crystallographic axis^^®. 

Perovskite compounds^^^ are in the chemical formula of ABO 3 , where A stands for alkaline 
or rare earth metal, B stands for transition metal (TM) and O stands for oxygen. Its 
structure is shown in Fig. 24(a) and (b). The ideal oxygen octahedral crystal field splits 
the TM d orbitals into Sg (with d^z'^_r 2 and dx 2 _y 2 orbitals) and t 2 g (with dy^, d^x and 
dxy) manifolds, separated by ~3eV in energy. Such crystal geometry usually is not good 
for band inversion. However, for a bilayer of the perovskite structure grown along [111] 
direction^^®, the TM ions are on a buckled honeycomb lattice as shown in Fig. 24(c). Based 
on tight-banding models and first-principles calculations, D. Xiao et al. have demonstrated 
that both t 2 g and Eg orbitals can generate topologically non-trivial states for TM ions on 
such structure. Especially, they showed that LaAuOs bilayers sandwiched by two LaAlOs 
substrates might possess a topologically non-trivial gap around 0.15 eV^^®. 

Various possible topological phases have been explored for the perovskite bilayer sys- 
tems^^®^^®°. Except the QSH phase, the QAH phase can be readily expected by breaking 
the TRS, as discussed before. Thanks to the presence of transition metal elements, different 
magnetic orderings may be realized in such systems with no need of additional doping. If 
the strong electron interaction among the d electrons is properly taken into account, we can 
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expect more interesting physics^^^ for such honeycomb lattices formed by bilayer perovskite 
structures. In a recent study, it was predicted^®^ by the LDA+Gutzwiller method^^®, that 
(111) LaCoOs bilayer grown on SrTiOs substrate may support the QAHE. Here, the inter¬ 
play between SOC and Coulomb interaction stabilizes a robust FM insulating state, which 
is a Chern insulator with non-zero Chern number. It was pointed out in this study, that 
the SOC splitting (which is usually not strong for 3d elements) can be greatly enhanced 
by the strong Coulomb interaction. It has also been revealed that the possibility of some 
intriguing phenomena, such as fractional quantum Hall effect, associated with the nearly 
flat topologically non-trivial bands can be found in Sg systems^®®. 

The perovskite bilayer structure can also be used to generate the spin-resolved QAHE 
state. Choosing elements adequately, one can fill Qg orbitals while leaving eg orbitals par¬ 
tially occupied, and the lODq gap (crystal-field splitting) will prevent t 2 g orbitals from 
participating in the low-energy physics. For an atomic sheet of buckled honeycomb lattice, 
the two B atoms in a unit cell contribute four Cg orbits in total. The band structure for these 
four orbitals is shown in Fig. 25(a) based on the tight-binding Hamiltonian with the hopping 
integrals given by the Slater-Koster formula^®^. Two dispersive bands are contributed from 
^)) which cross each other linearly and form the desired Dirac cones at K and K' 
due to the honeycomb structure. 

A candidate perovskite material which generates d® electrons is LaAuOs, wherein Au 
takes the -|-3 charge state. In order to generate an AFM exchange field, one can sandwich 
the atomic sheet of buckled honeycomb lattice with d® electrons by two thick layers of 
insulating perovskite material grown in the [111] direction. A candidate material for the 
substrates is LaCrOs, a Mott insulator with C-type AFM order, which induces opposite 
exchange fields on the two sublattices in the buckled honeycomb lattice. The large Mott 
gap of LaCrOs prevents undesired influences on the Dirac behavior of d®’® electrons except 
for the exchange field. Since La is chosen as the common A-site cation, the lattice mismatch 
between substrate and the targeting atomic sheet is minimal. First-principles calculations 
were performed in terms of VASP^^^ code with a supercell of (La 2 Au 206 )/(La 2 Cr 206)5 shown 
in Fig. 25(b). As displayed in Fig. 25(c), the black bands are mainly from the substrate 
LaCrOs with a large gap of 3eV, characteristic of a Mott insulator. The bands inside the 
Mott gap are contributed by the Cg electrons of Au+^ ions. There are two almost fiat bands 
and two dispersive bands, as captured by the tight-binding Hamiltonian (see Fig. 25(a)). 


68 


(a) (b) 



FIG. 25: (Color online) (a) Noninteracting tight-binding band structure of the eg system, 
where each band is labeled with the corresponding electron conhguration if it is occupied. 
Inset shows the two eg orbits, (b) Schematic of a unit cell (La2Au206)/(La2Cr206)5. (c) 
Band structures based on first-principles calculation. Figures from Ref. 161. 

Gaps are opened at K and K' since the Dirac electrons now feel the AFM exchange field from 
LaCrOs and SOC. The two gaps are different in magnitude and given explicitly by 2(M-|-A) 
and 2(M-A) respectively (see Fig. 23(b) for V = 0), from which one obtains M = 0.166eV 
and A = 32.91meV for this system. SOC effect is quite large in the present case, due to the 
orbital hybridization between the two eg orbits in the new frame along the [111] crystalline 
direction. All the states in Fig. 25(b) are doubly degenerate due to the product symmetry 
of TRS and IS. If a gate electric field is applied along the [111] direction, the Au ions in the 
buckled honeycomb lattice should experience staggered electric potentials. By increasing the 
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electric potential, the system should be driven into the spin-resolved QAHE, as summarized 
in Fig. 23(b). The required electric field is on the order of O.IV/A. As the maximal non¬ 
trivial gap is given by 2A ~ 600K, this topological state in the LaCrOs/LaAuOs superlattice 
may remain stable even around room temperature. Interaction driven QAHE phase can also 
be found in SrTiOs/LaCoOa/SrTiOa quantum well as discussed in detail in Ref. 218. 

It is worth noticing that, although the edge current of the spin-resolved QAHE state is 
fully spin polarized, the total system does not show any net spin magnetization, since the 
numbers of occupied states in the two spin channels of the buckled honeycomb lattice are the 
same, and the substrates are AFM insulators. The spin polarization of the edge current can 
be reversed by the gating electric field, which is superior for future spintronic applications. 


V. DISCUSSIONS AND FUTURE PROSPECTS 

At this stage, the quantum Hall trio (namely, quantum Hall effect, quantum spin Hall ef¬ 
fect, and quantum anomalous Hall effect) has been completed.Looking back to history, it 
was a long term exploration towards the final realization of QAHE. First, the understanding 
of the Berry phase mechanism of intrinsic AHE was an important step, because it established 
the link between the AHE and the topological properties of electronic structure. Second, 
the rapid progresses in the field of topological insulators strongly stimulated the studies on 
QAHE. Third, first-principles calculations play important roles in predicting and selecting 
candidate materials and systems. This also establishes a nice example that theoretical and 
computational studies can in some cases really go ahead of and greatly help experiments. 
Finally, the well-controlled growth and magnetic doping of high quality samples and trans¬ 
port measurement at ultra-low temperature make the experiments possible, thanks to the 
progress of modern experimental technology. All these aspects combined together finally 
gave birth to the QAHE. 

On the other hand, the successful realization of QAHE is not the end of the story. Instead, 
we believe this success will inspire more extensive research in this field. Two issues become 
immediately important for future studies: 

(1) How can we realize a higher plateau with a Chern number larger than 1? 

(2) How can we increase the temperature range of QAHE (it is now observed only in the 
tens mK range)? 
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From the theoretical point of view, the QAH states with higher Chern numbers can be 
derived in a way similar to that described in section VLB if the magnetic exchange field can 
be made largerLarger exchange field will, in general, generate more band inversions 
among sub-bands and thus induce higher Chern numbers. While in practice this idea is very 
difficult, because the system with large enough exchange field will usually become completely 
metallic before the second plateau state appears. An alternative way to realize the QAHE 
with a high Chern number is to start with materials that have multiple band inversion 
regimes in the BZ. An interesting proposal is the magnetically doped topological crystalline 
insulators (Pb,Sn)(Te,Se), which have band inversions around all the four T points in the 
BZ^^h Similar to the situation in Bi 2 Se 3 family compounds, here the band inversions appear 
between p-bands with opposite parity, and the Van Vleck mechanism also works, which can 
lead to a stable FM semiconductor ground state after magnetic doping. Therefore, for thin 
films of this system with proper thickness, it is possible to realize the QAHE with a Chern 
number between ±4^^^. 

The experimental observation of QAHE in thin films of Cr-doped (Bi,Sb) 2 Te 3 requires an 
extremely low temperature (approximately 30 mK), which hinders realistic applications^^^. 
This is presumably due to the small band gap opened by exchange coupling and the low 
carrier mobility of 760 cm^/Vs in the sample. Therefore, searching for realistic systems with 
non-trivial band structures, strong exchange coupling and high sample quality is essential 
to realize the QAH effect at a higher temperature. Liu et al. proposed to realize QAHE 
in a more conventional diluted magnetic semiconductor, the Mn-doped InAs/GaSb type H 
quantum well, which has a higher carrier mobility of 6000 cm^/Vs and a predicted FM 
transition temperature of 30 Other theoretical proposes for high temperature QAHE 

systems are based on the stanene and germanene materials. These two materials are re¬ 
cently predicted to be QSHI with a large insulating gap^^^. The QSH phase appears due 
to an energy inversion between the Sn(Ge)-s and Sn(Ge)-p 2 ; 2 / bands, which have opposite 
parities and lead to a sizable bulk gap of 0.3 eV. In Ref. 223, S. G. Wu, et al, proposed 
that half-passivated stanene and germanene become Ghern insulators with energy gaps of 
0.34 and 0.06 eV, respectively. The proposed EM order in this system comes from the 
unpassivated sublattice that exhibits dangling bonds. Strong coupling between the spin- 
polarized dangling bond states and the inherent inverted bands opens a considerable energy 
gap. The estimated mean-field Gurie temperature is as high as 243 K and 509 K for Sn 
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and Ge lattices, respectively. In spite of those theoretical proposals, it is still challenging to 
reach QAHE at higher temperature experimentally. To achieve this goal, searching for new 
host materials with simpler electronic structure, easier synthesis and less extreme growth 
conditions is necessary. 

The QAHE and Chern insulators can be also helpful for the possible realization of 
Majorana fermions (ME)^^^, which are exotic charge-neutral particles that are their own 
antiparticles. In recent years, Majorana fermions have been intensively studied in the 
held of condensed matter physics due to their importance to both fundamental science 
and potential applications in topological quantum computing. Theoretically, Majorana 
fermions have been predicted to exist in several spin orbit coupled superconducting systems, 
such as superconductor-topological insulator interfaces'^, semiconductor-superconductor 
heterostructures^^®^^^^, nanotube-superconductor devices^^®, half-metal-superconductor in¬ 
terface^®, etc. Experimentally, great efforts have been made in this direction, and some 
important evidence has been successfully observed^^®’^^®. On the other hand, in Ref. 75, a 
different way to generate Majorana fermions is proposed in a hybrid device that consists of 
a QH or QAH insulator layer in proximity to a fully gapped superconductor layer on top. In 
this case, it can generally be expected that a chiral topological superconductor phase with 
an odd number of chiral Majorana edge modes will exist. Compared to the QH system in 
the proposal, the QAH system has the advantage that it does not require external magnetic 
field, which makes the realization of superconductor proximity effect much easier. The Ma¬ 
jorana zero mode in a model combined with a QAHE system and an s-wave superconductor 
is also discussed in Ref. 231. 

In summary of the present paper, we have reviewed the history of achieving the QAHE. 
In 1880, E. H. Hall found the anomalous Hall effect, 17 years before the discovery of the 
electron. In the following 100 years, however, the developments in this field, particularly for 
the understanding of AHE, were slow. There were debates on whether intrinsic or extrinsic 
contribution is dominant. The breakthrough was the discovery of the QHE in 1980, from 
which people learned that there exist new quantum states characterized by the topology of 
electronic structures. The concept of the Berry phase, established in 1984 by M. V. Berry, is 
the solid base of various topological quantum states. We have reviewed how the Berry phase 
is related to the AHE, QHE, QAHE, QSHE, 3D TIs and 3D Weyl/Dirac semimetals. Such 
a unified view has led us to a unified scheme, called the Wilson loop method, in calculating 
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various topological invariants, including the Chern number, Z 2 index, mirror Chern number 
and Fermi surface Chern number. This method is efficient in rigorously determining the 
topology of real materials. However, in practical calculations, the band inversion mecha¬ 
nism is much more intuitive and efficient in predicting, selecting and designing candidate 
topological materials. Finally, we have listed the necessary and complete four conditions 
for realizing the QAHE. Following these conditions, we have demonstrated several ways to 
achieve the QAHE. The theoretical prediction for the QAHE in magnetic-ion Cr and Fe 
doped Bi 2 Se 3 and Bi 2 Te 3 family TIs has been recently confirmed by experimental observa¬ 
tion. The key difference of this scheme from other proposals is the FM ordering of Cr and 
Fe dopants induced by the van Vleck paramagnetism, which makes this proposal successful. 
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